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Abstract 

We investigate quasi-hermitian quantum mechanics in phase space using standard deformation quan- 
tization methods: Groenewold star products and Wigner transforms. We focus on imaginary Liouville 
theory as a representative example where exact results are easily obtained. We emphasize spatially 
periodic solutions, compute various distribution functions and phase-space metrics, and explore the re- 
lationships between them. 
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1 Introduction 



Superficially non-hermitian Hamiltonian quantum systems are of considerable current interest, especially in 
the context of PT symmetric models [3j [13] . For such systems the Hilbert space structure is at first sight 
very different than that for hermitian Hamiltonian systems inasmuch as the dual wave functions are not just 
the complex conjugates of the wave functions, or cquivalently, the Hilbert space metric is not the usual one. 
While it is possible to keep most of the compact Dirac notation in analyzing such systems (see Appendix 
E), in the main body of this paper we will work with explicit functions and avoid abstract notation. Our 
goal is to expose the underlying mechanisms (as in [5J [3]) rather than to hide them. 

Our discussion is focussed on a system with potential exp (2ix). This model, as well as its field theory 
extension, is of interest for applications to table-top physical systems (2J [4] and to deeper problems in string 
theory [l9j|20j. We will not discuss those applications here, but rather we will simply develop the phase- 
space formalism for the point particle model. Wc believe this formalism will be helpful in understanding 
the applications cited, as well as others. Other recent work along these same lines can be found in [161 117j . 

2 Imaginary Liouville quantum mechanics 

Consider "imaginary" or "periodic" Liouville quantum mechanics as governed by the apparently non-hermitian 
Hamiltonian 

H = p 2 + m 2 e 2lx (1) 

Without essential loss of generality, we take to = 1 in most of the following. Obviously, this is a "PT 
symmetric" model. But more to the point, this is actually a "quasi-hermitian" theory [18) with a real 
energy spectrum, as explained in [8J and as we shall clarify further here. We will analyze this system in 
phase-space using the methods of deformation quantization. 

2.1 Eigenfunctions 

But first, let us briefly review the position representation Schrodinger eigenvalue problem for this system 
(see e.g. [5]). With x on the real line and with the condition that the wave functions remain bounded, the 
corresponding Schrodinger equation has energy eigenvalues given by all real E > 0. The eigenfunctions are 
just Bessel functions, </ ±v /£ (e IX ). These are doubly degenerate when \[~E ^ n € N, but they merge into 
a single, nondegenerate eigenfunction when E — n 2 . For these nondegencrate cases the eigenfunctions are 
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27r-pcriodic in x, and the solutions are the analytic Bessel functions J n 

/ d 2 



n4 ^ = i"tf +e j Jn(eB) (2) 
-i 00 ('_1"l fc „2ifc:E 

J n ( e ix ) = —e inx Y , , \ ' . , , n = 0,1,2,- •■ (3) 
v ; 2™ ^ fc! (fc + n)! 4 fc ' ' ' w 

We will consider only such periodic solutions here, since in our opinion this is the most interesting situation. 
We note that these periodic eigenfunctions are superpositions of only right-moving plane waves. In fact, 
in this periodic situation the discrete energy spectrum is precisely the same as would be found for particles 
moving freely on a circle but restricted to non-negative momentum. 

2.2 Dual polynomials 

The periodic Bessel functions and their complex conjugates do not form an orthonormal set on the circle. 
To obtain an orthonormal set of functions it is necessary to combine { J n (z)} with an associated set of 
polynomials in z , {A n (z)}, the so-called Neumann polynomials. These are dual to {J n (z)} on any 
contour enclosing the origin z = 0. in the following sense: Ai ( z ) Jk ( z ) = $j,k- So on the circle, we 

have 

dxA j (e i *)J k (e ix )=8 j , k (4) 
'e given explicitly bjQ 



L !^ J (r, _ h- IV p 2ikx 

A (e ix ) = 1 , At (e ix ) = 2e~ ix , A n (e ix ) = 2 n ne~ mx ^ * '-^ (5) 

k=0 

l n (^) = (-^ + ^-n^A n (e-) (6) 

As indicated, the {^4 n } obey inhomogeneous modifications of Bessel's equation, where the inhomogeneity is 
cither X n (z) oc z 2 for even n or T n (z) oc z for odd n, according to 

e n e 2lx for even n > 
In (OH " (7) 

2ne lx for odd n > 

While © is inhomogeneous, nevertheless the usual proof of orthogonality between pairs of non-degenerate H 
eigenfunctions and their duals goes through because the inhomogencitics arc orthogonal to the eigenfunctions. 

dx Xj (e ix ) J k (e ix ) = , for all j, k G N (8) 

o 



1 For convenience we have modified the usual notation of the associated polynomials as given in [T] 1151 1221 , namely O n , and 
have defined A n (z) = e n z O n (z) where eo = 1 and e n = 2 for n 0. 
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3 Phase space distributions 

We next compute Wigner transforms of various function bilincars for the 27r-periodic Bessel/Ncumann sys- 
tem. For a system so-defined, on a circle, momentum is quantized. Thus we would expect that the (x,p) 
phase space is not the usual R 2 nor even S 1 x R, but rather that it is reduced to S 1 x Z. In fact, the 
periodic energy eigenfunctions of the imaginary Liouville Hamiltonian consist of superpositions of positive 
momentum plane waves, so we would also expect not to need the Z<o momentum sector at all. Well, both 
expectations are almost true. But not quite. We shall see below to what extent these expectations are 
born out. 

3.1 Eigenfunction WFs 

As a first step, we remind the reader about the structure of real (diagonal) Wigner functions (WFs) made 

from 27r-periodic plane waves. They are just Kronccker deltas, with p £ Z as expected. 

1 f*« 

e„ (x,p) = ^J o 4>n (x - y) K (x + y) e 2lyp dy (9) 

= 5 n>P (10) 

upon choosing <p n (x) = exp (inx) . By analogy, for the Liouville eigenfunctions ?/>„ the WFs are again 
manifestly real, and again have support for p G Z as expected, as given by 

f n (x,p)=±-J ^ n {x-y)i^{x + y)e 2 ^dy (11) 

f_-\\P~ n P ~ n e 2ix(n-p+2k) 

= 4~p ^ fc! (n + k)\ {p -k)\(p-k- njl ^ 
where the sum results from taking tp n (x) = J n (e ia: ) as given by the series in ([3]). Note for the Liouville 
case, as opposed to the free particle case, the support in p is infinite: This particular /„ (x,p) is non-zero 
for all p > n > 0. Another way to write these WFs makes use of the associated Lcgcndrc functions! 2 ] 

(-1)" Asin2a;\ p 

/„ (x,p) = — — - — LegcndrcP (p, -n, i cot 2x) (13) 



p\ (p — n)\ \ 2 

In particular, for any point in the reduced phase space, (x,p) 6 S 1 x Z, we have 

fo (x,p) = 1 x S Pj0 - ^ (cos 2a;) x S P:1 + i (2 + cos Ax) x 8 pa H (14a) 

/i (x,p) = — x <J p ,i - i (cos 2a;) x 6 P:2 + ^ (3 + 2cos4a;) x 5 Pj3 H (14b) 

fs(x,p) = ^ x 6p, 2 - (cos 2a;) x S Pj3 + (4 + 3 cos 4a;) x 8 pA + ■ ■ ■ (14c) 

h (x,p) = — x $p,3 - 2TT32- (cos2x) x S pA + (5 + 4 cos 4a;) x 5 P , 5 + ■ ■ ■ (14d) 



A form which facilitates a continuation to non-integer p, should anyone wish to do that. 



But now we have more functions at our disposal, namely the Neumann polynomials, so we may build another 
set of WFs for comparison to those in (JTUJ) and (fT2")l . 



3.2 Dual WFs 

For dual functions Xn the WFs are also manifestly real, with p G Z, as given bjf| 

f n (x,p) = —J Xn (x - y) Xn (x + y) e 2lyp dy (15) 

|n/2j |n/2j 

= 4%2 ^ ^ ( n-fc- 1) K n- - 1) ! e2M; - fc) (M) 

fe=0 i=0 

where the sums result from taking Xn (x) = A n (e IX ) as given by the series in ([5]). Note the support in p is 
now finite: This particular /„ (x,p) is non-zero for < n — 2 [^/2J < p < n.. That is to say, < p < n for 
n even, and 1 < p < n for n odd. 

Resolving the constraint set by the Kronccker delta in (|16p eliminates one sum and further restricts the 
range of the other to yielqf 

min(|n/2J,n-p) ,,,,, ^ 

£(*,p) = 4?n» £ ( ""fc,7„l ( _ fc 7 ^ "^-^ ( 17 ) 

fc=max(0,n-p-Ln/2j) " 

In particular, again for any phase-space point (x,p) £ S 1 x Z, 

/o = <5p, (18a) 

/ 1 (x,p)=4V (18b) 

/ 2 (x,p) = 4 x (5 P;0 + 32 (cos2x) x <5 p ,i + 64 x <5 P:2 (18c) 

/ 3 = 36 x dp.j + 576 (cos2x) x 5 P , 2 + 2304 x 6 P>3 (18d) 

Thus, at any given momentum level, we find the same set of functions of x (i.e. cos(2kx)) no matter whether 
we consider {/„} or j/nj- 

There is a basic orthogonality relation for WFs and dual WFs. 

Tj- £ fk(x,p) fn(x,p) = S k ,n (19) 

This follows in a straightforward way from the bilinear structure of the Wigner transform and from the 
orthogonality of the wave functions and their duals. There is also a corresponding pseudo-local relation on 



3 To dispel any confusion about our conventions for f n , the complex conjugation in H15H is different from that in (J9j and mil l 

just because we chose in [8] [9] to define the duals such that <5 n j. = i \n (&) rfk ( x ) dx without any explicit conjugations. 
4 The limits on the sum may also be written as min ( [™/2j , n — p) and n — p — min ( L«/2J , n — p). 
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the phase space that involves the Groenewold star product (see (|61|) and (f6"5"]) below). The form of these 
results can be seen most easily through the use of formal density operator methods, as in Appendix E. 

Perhaps it is useful to present the specific examples of /„ and f n , for n = 0, 1, 2, 3 and for p = 0, 1, 2, 3, in 
Table form. This facilitates checking (fill]) for these few cases, and illuminates the orthogonality mechanism. 







WFs & Dual WFs (non-zero 


values) 




p = 


p=l 


p = 2 


p = 3 


fo, fo 


1 , 1 


— ^ cos 2x , 




1152 C ° S6X 128 C ° S2X '° 


fljl 




I* 


— — cos 2x , 
16 


1 „ 1 n 

— cos4x+— ,0 


h,f2 


, 4 


, 32 cos 2x 


64' 64 


-^cos2x,0 


f3, $3 




, 36 


, 576 cos 2x 


1 , 2304 
2304 ' 



4 Wigner transform of the bilocal metric 



It is explained in [8] - as well as in the classic literature on the subject - how a scalar product for a 
biorthogonal system such as {Ak,J n } can always be written as an integral over a doubled configuration 
space involving a "bilocal metric" K (x,y). 

(0, V) - J J (x) K (x, y) i> (y) dxdy (20) 
4.1 Bilocal <-> phase space 

When a scalar product is so expressed as a bilocal bilinear form then it is naturally and very easily re- 
expressed in phase space (which we suppose to be M. 2 in this paragraph) through the use of a Wigner 
transform [IT] . 

f 4>4 , ( x , p ) = ±-Jtp(x- \y\ <j>(x+ \y^j e*w dy (21) 
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We have chosen the normalization here so that for p on the real line 

poo 

ip(x)</>(x) = f^, t> {x 1 p)dp (22) 



More generally, Fourier inverting (|2ip gives the point-split product 

</> (x) i> (y) = J°° jto-'bUt (^-A dp (23) 
Thus the scalar product (|20|) can be re-written as 

0, = llR (x,p) U<t> (x, p) dxdp (24) 



where the phase-space metric is the Wigner transform of the bilocal metric. 



and inversely 



K(x,y) = ^J°° e**-y>R (^>p) dp 
In a more abstract notation (as in Appendix E) the form of ([24)1 is [4>, ip) = Tr (R (0|) = Tr I \if; 



(26) 



4.2 Liouville dual metric 

The preceding results are quite general. To be more specific, for 27r-periodic dual functions of imaginary 
Liouville quantum mechanics, the scalar product was shown in |8] to be , 2 \a ff Q ™ dxdyxj (x)J (x, y) Xk (y) 
= 5j.k where 

oo 

J (x, y) = Jo (e~ ix - e*) = J (e~ ix ) J (e*») + 2 ^ J n [e~ ix ) J n (e*) (27) 

n=l 

Again, just a Bessel function. Or, re-expressed in a form which is immediately useful in the following, 



J{x,y) = J Q (-2ie^- x ^sm(^±^ 



(28) 



Up to a normalization the corresponding metric in phase space is given by the Wigner transform of this 
bilocal| 

A bit of care is needed since the Wigner transforms f n (x,p) on which this metric will act are actually 
defined so that a dual function \ pl a Y s the role of <f> and a conjugate dual function \ plays the role of tp in 
the above (compare (f2Tj) to (fl~5|) ). Thus, acting on /„ (x,p) the metric would be the Wigner transform of J 
as above, only with first and second arguments interchanged. We also adjust the normalization here (and 



5 To take the free particle limit, the parameter m in JTJ must first be restored. See Appendix A. 
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again later, in (|33p ) to take into account our conventions and the fact that we are dealing with 27r-periodic 
functions (see Appendix C). In view of all this, we finally obtain a dual phase-space metric given by 



i r 2 * , . o,„..„ i r2n 



R(x,p) = —J J(x + Wl x-w)e 2iwp dw = — J J Q (~2ie~ lw ainx)e 2iwp dw 



Hence the simple final answer. 



I ■ 2 \v 

i? (x,p) = j— for integer p > 0, but vanishes for integer p < (30) 

! ) 

An equivalent operator expression can be obtained by the method of Weyl transforms. (See Appendix F.) 
Another way to characterize R(x,p) is to note that it satisfies the differential-difference equation 

pd x R (x,p) = sin (2x) R (x,p — 1) (31) 

even when p = 0, since d x R (x,p = 0) = as well as R(x,p = —1) = 0. (This should be compared to p41[) 
given below. Note that R actually corresponds to in that later discussion.) 

In fact, there is an obvious continuation of (|30[) to all real p, or even to complex x and p. Namely 



~ , (sin 2 *)* , , 

with its manifest zeroes and singularities (poles and cuts). This continuation also transparently satisfies 

(ED. 

Taking into account all our conventions, we may now express the normalizations of pure states in terms 
of the dual WFs and the dual phase-space metric as 

i r ~ ~ 1 n f 2n (sin 2 x) P ~ 

£n = ^±R(x,p)f n (x,p) = —J2 V - ' fn(x, P )dx (33) 

where as usual £q = 1 and e„ = 2 for n > 0. It is tedious but straightforward to use (JT7J) to check this 
normalization and confirm that the dual metric does its job. More importantly, (|33[) is consistent with (|19[) 
for the simple reason that 

oo 

R(x,p)=J2 e kfk(x,p) (34) 
This in turn follows from the expansion of the bilocal metric in terms of Bessel bilinears, in (|27D . 



5 Homogeneous versus inhomogeneous ^-genvalue equations 

Were the dual functions just the complex conjugates of the wave functions, the two types of WFs that we 
have defined would be identical, /„ (x,p) = f„ {x,p)\^ =x , but this is obviously not true for the case at hand. 
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5.1 WFs as eigenfunctions 

As an important distinguishing feature, for the first of these WFs we have 

H*f n = n 2 f n = f n *H (35) 

where the associative Groenewold star product operation [T2] is (h = 1) 

★ = exp d x dp - ^ dp d x ^j (36) 

Thus the wave function WFs are energy -fcgenfunctions. 

The -^genvalue equation (|33|) for the WFs is not only a distinguishing feature but can actually be used 
to define the /„ as real functions, without any prior knowledge of the underlying wave functions. Taking 
imaginary and real parts of (|35|) we obtain an equation like (|3~lj) , namely 

pd x fn (x,p) = sin (2x) f n (x,p- 1) (37) 

as well as 

p 2 - jd^j /„ (x,p) + cos (2x) f n {x,p - 1) = n 2 f n (x,p) (38) 
Using (|3"T|) twice, the second derivative in ([55]) becomes 

^9 » / x 2cos(2x) „ , _ sin 2 (2a;) , . „. ,„„. 

dlh (x,p) = K -^f n (x,p- 1 + v /„ (x,p- 2 39 

p pCp-1) 

Hence (J38J) reduces to a second-order difference equation in the momentum. 

(p 2 - n 2 ) f n (x,p) + cos (2s) /„ (s,p - 1) - - 2) = (40) 



iitiony 



We may solve this second-order equation by forward recursion under the condition^ that /„ (x,p < 0) = 0. 
The resulting WFs have support only for non-negative integer p. We find that /„ (x,p < n) = 0, f n (x,p = n) 
is arbitrary, and all f n (x,p > n) are uniquely determined by (|40|) in terms of our choice for /„ (x,p = n). 
For example, for the ground state n = 0, the choice /o (x, 0) = 1 immediately reproduces the terms in (|14ap . 
Similarly the choice (with no x dependence) 

fn(x,p = n) = (41) 
4™ (n!) 

reproduces the series (|12[) . The choice for /„ (x,p = n) must be independent of x so that /„ (x,p < n) = 
is consistent with (|37p . Thus the WFs are determined by the -^rgcnvalue equation. 



J Allowing support for negative integer p leads to singularities at p = and p = 1. Therefore we rule out this possibility. 
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5.2 Entwining the dual metric 

It follows immediately from (|34[) and (|3"5)) that 

H*R(x,p) = R{x,p)*H (42) 

which amounts to just (|3Tj) . (This should be compared to (|135[) given below.) The equivalent operator 
statement is obvious, and follows directly from the Weyl correspondence. (Sec Appendix F.) As direct 
verification of (|4*2"j) . we compute 

H*R(x,p) = — R(x,p) = R{x,p)-kH (43) 
2 sin x 

The relation (|42"]) is actually a special case of the "two star equation" given in [7], Eqn(73). In the language 
of that paper, (|42[) is the ultra- local version for which T (x,p; X, P) — S {x — X) 8 (p — P) R (x,p) and one 
Hamiltonian function (H) is the complex conjugate of the other (TC). 

5.3 Inhomogeneities for dual WFs 

For the dual WFs (fTS")) . by direct calculation using standard methods (as in an Appendix B) we find 
In * H = M p + l -id x ^j + e 2lx e a ^j i- J 2 J x^(x-y) X n (* + y) e 2i ™dy 

= ^ jj e 2 ^ X^(x-y) (~ (d y + d x ) 2 + e^+y^j Xn (x + y) dy 
Then from (|6j) we have 

2ne l( - x+v ^ for n E N odd 

Z*H = n 2 Z+±- f * dy e 2 ^x^(x-y)< 
27T Jo 

Taking Xn (%) — A„ (e M ) and using ©, this gives 

fro 4fefc! 2 ^o | eneMM 

= ^ + 2«, V - *"v» - { 2 " e ^~ 

i 1 4"fc! c 2ia;A 

fe=0 £ n 6 £>2p+2/c,n-2 



(44) 



(45) 



,2i(x+ a ) for n eleven 



(46) 



That is to say, in complete analogy with the dual wave functions, the dual WFs obey inhomogeneous 
•^-genvalue equations. 

In particular, for the dual ground state WF 

f *H = e 2ix 5 p ,- 1 , H*f = e- 2iX 5 p ^ 1 (47) 
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and more generally 



f n *H = n 2 f n + < 



2 (L»/2J+p)! 4 l+p e 2i(X+p) x 
(|n/2J-p)l 



for odd n, if [n/2\ > p > 0, otherwise 



2n ( [„/2j -p-i) ! 4 1+ Pe 2i ^ 2+ ^ X for even n > 0, if L"/2J > 1 + P > 0, otherwise 



(48) 



2 (L"/2j+p)! A l+p -2»(l+p)a 
" (Ln/2j-p)! 4 e 



for odd n, if ^tt, / 2 J > p > 0, otherwise 



#* /n = n?fn + < 



(49) 



2n ( Ln/2j -p-i)! 4l+Pe ~ 2 ' (2+p)a: for cvcn " > 0, if L«/2J > 1 +P > 0, otherwise 
As an alternative derivation, this can be checked by acting directly on the dual WFs as explicit sums, (fl7|) . 

These equations for the dual WFs do not determine them uniquely, when they are only required to be finite 
real solutions to the implied differential-difference equations (again with /„ [x,p = n) independent of x) since 
we could always add to any /„ a solution to the homogeneous equation, namely an /„. However, if we require 
the stronger conditions that we seek solutions with finite momentum support, such that /„ (x,p < 0) = 
and /„ (x,p > n) = 0, then the solutions arc uniquely determined, up to normalization. In fact, the form of 
the inhomogeneities is also fixed by these requirements. This is similar to what happens in the analysis of 
the dual wave functions, which may be constructed without knowledge of the wave functions if we require 
that they be finite polynomials in cxp (— ix). In that analysis the form of the inhomogeneities is also fixed 
(see [8] and [9]). 

In (|47[) . and in (|4"5)) and (|4T)|) for other even n, we see that the inhomogeneity has support for p = —1, 
even though /„ has support only for non-negative p. The star product with the Hamiltonian has spread out 
the distribution, in a typical quantum fashion, just slightly into the realm of negative p. Thus our original 
expectation that we should be able to ignore all integer p < was not quite correct. 



5.4 (H) from dual WFs 

For imaginary Liouvillc QM the difference between (TJ5|) and (fJ5|) is 

2 ^ 2 ([»/2j-p)! 4l+P sin ( 2x U + P)) for odd n > if L«/2J > P > 

f n *H -H-kf n = < 



(50) 



Ain an%i-t-iy. Al+P sin ( 2x ( 2 + p)) for cvcn 71 > °' if L«/2J > I + P > o 

The peculiar structure of the RHS combines with that of the metric (|30[) to establish specifically for the 
Liouville model that the expectation of the Hamiltonian is real within the phase-space framework. 

I 



= 



2tt 



R(x,p)(f n *H -H*f r 



(51) 
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This holds just because = (sin 2 x) P sin (2x (1 + p)) dx as well as = J Q ™ (sin 2 x) P sin (2x (2 + p)) dx, 
for integer p > 0, and because the vanishing of the metric for negative p conveniently eliminates the p = — 1 
possibility for even n. 
Indeed, for all n 



n 2 



j{ R(x,p)(f n (x,p)*H) Y[R(x,p)(H 

J x,p J x,p 



fn (x,p) 



R{X,P) fn (X,p) M R{x,p)f n (x,p) 

x,p Jx,p 



(52) 



despite the inhomogeneities in (|48[) and (|49|) . because = f ™ (sin 2 x) P exp (2ix (1 +p))c?x as well as = 
Jq 77 (sin 2 x) P exp (2ix (2 + p)) dx for integer p > 0. Again the vanishing of the metric for negative p eliminates 
the possibility of anguish at p = — 1 Q 

However it must be said that (|51[) can also be established more generally by the method of combining 
(|42[) with the "Lone Star Lemma" |23j which allows us the option of inserting or removing a single * from 
the phase space summand/integrand, hence to write 

^-£r(x, P ) (Z*H) =^£R(x,p)*U*H=^£ i H*R(x,p)*fZ 

= ^-'£R(x,p)*^*U = ^-'£R(x,p) (h*Z) (53) 

6 A sesquilinear star product and bracket 



There is a compelling formal way to express the entwining relation (I42p between the dual metric and the 
Hamiltonian which suggests it would be appropriate to abbreviate quasi-hermitian to q-hermitian in the title 
of this paper. Certainly, the structure of (|42p brings to mind the deformations of commutators and braiding 
relations such as occur in q- algebras and quantum groups [6]. To pursue that statement, let us rewrite the 
entwining relation as 

H*R(x,p) = R(x,p)*KHK (54) 

where K is the anti-linear operation of complex conjugation. Since R is in fact a real function, this may 
also be written as 

H*KR(x,p) = R(x,p)*KH (55) 
By defining a "sesquilinear star product" 

© = * K (56) 

7 Otherwise the RHS of, say, I@7J|, would be troublesome when multiplied by (sinx) 2p and integrated over x, although 



l P =-l 



not impossible to handle (cf. Cauchy's principal value prescription). 
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and a corresponding modification of the Moyal bracket [14] , the entwining relation becomes simpljo 

= \H, r] = H®R - R®H (57) 
L J© 

A form of sesquilinearity is evident, although left<->right ordering dependent, in the properties 

C® {aA + f3B) = (C®A) a + (C®B) f3 = a (C®A) + ff (C®B) (58) 
(a A + (3B) ®C = a (A®C) + (3 (B®C) = (A®C) a + (B®C) (59) 

for constants a and f3. Also note that -kK ^ K*, rather 

★ = IK * K = exp ( ~ d x dp + | d p d x J (60) 



This "star-bar" is again an associative but non-commutative product, like (|36[) . and just amounts to flipping 
the sign of the deformation parameter, h (= 1 throughout this paper). 

7 Other phase space distributions 

We may also take star products of the various WFs. This leads in a routine way to a larger class of 
distributions on the phase space, although the details for the present circumstances are rather novel not only 
because the functions are periodic in x but also because \n 7^ ipn- This means that in general there is a 
third class of "hybrid" WFs involving mixed bilinears in the wave functions and their duals. 

7.1 Hybrid WFs 

By direct calculation, with p £ Z, we find 

fk(x,p)* f n (x,p) = g n (x,p) 6k, n (61) 



(x,p) = ^ jj {x - V) Xn (x + y) e 2 ^dy (62) 



* (x, P ) = A S~T\, e Mp - n)X t! J»r*~Jl { 7? k ^ X (63) 



4P-" (n - 1)! 4 2fe fc! (p + k)\(p + k- n)\ 



8 For the less symmetric Voros product <l = exp yi d x 9 p j (see the Appendix in |21| ) an analogous construction is not quite 
so elegant, since in that case the metric would not be a real function. 



15 



the latter from the series and ([5]). So g n (x,p) 7^ if and only if p > n — \n/2\, and thus has infinite 
support in p similar to /„. In particular, again for any phase-space point (x,p) £ S 1 x Z, 

go(x lP ) = ^P^ e 2lpx for p>0 (64a) 
4p (pi) 

g2 = 4P (p lm P D! ( 16 ( p2 - e " (p ' 2)3: - e2 " p ) for (64c) 

etc. Without further calculation we also find 



f k (x,p) *fn (x,p) = g n (x,p) d k , n (65) 



9n(x,p) = — \n (x - y) ip„ (x + y) e 2lyp dy (66) 

Z7T 



II 



as consequences of the previous relations, upon making use of /„ * f n = /„ * f n . Summing and integrating 
over p and x the RHS of either ([61]) or (|65|) gives back (fT9|) . A sum over integer p suffices in this case|£ 
These hybrids also obey a form of the completeness relation that follows from the corresponding statement 
for wave functions and their duals. 

oo 

=Y, A n M Jn (Z) (67) 

However, the Wigner transform requires some regulation to give a precise meaning to X^fcLo 9 k ( X 'P)- 
Finally, we find the remaining simple star products 

g k (x, p) * g n (x, p) = g n (x, p) S k>n (68) 

9k{x,p)~k f n (x,p) = f n (x,p) Sk,n (69) 
fk(x,p)*g n (x,p) = f n (x,p) S k ,n (70) 

The first of these is of the usual esthetic form that holds for non-hybrid WFs in familiar hermitian quantum 
mechanical systems (e.g. the free particle WFs satisfy e k (x,p) *e„ (x,p) = e n (x,p) S k . n )- However, unlike 
in that non-hybrid hermitian situation, the hybrid g n is not real. A direct consequence of the second of 
these star products is that all the g n are right-*- mapped by the dual metric into /„. 

g n (x,p)*R(x,p)=e n f n (x,p) (71) 



9 However, in other situations it may be necessary to sum over semi-integer p, i.e. all p such that 2p £ Z. While e n , f n , f n , 
and g n only have support for integer p, there are "non-diagonal" variants of these distributions, namely e k n , fk.n, fkny an d 
gk,n f°r k 7^ n (see Appendix D), which have support for semi-integer p when k — n is an odd integer. In dealing with such 
nondiagonal cases, some care is required to respect the positivity of the metric. 



16 



This follows immediately from 

Once again, the form of the various star products of WFs can most easily be seen through the use of 
density operator methods, as in Appendix E. However, for devotees of phase-space methods, we make the 
following technical points. The results are obtained directly from the structure of the Wigner transforms 
(fTTjl . (fT5|) . and (p)2"|) , and the orthogonality of the wave functions and their duals. For example, consider 
(|6T|) . Applying the star product to each of the integral representations (fTTj) and (fT5|) . and carrying out the 
requisite variable shifts on the integrands, the resulting coupled integrals J 2w dyi dyi ■ ■ ■ can be split 
into uncoupled integrals d (j/i + 2/2) • • • x \ d (yi — 2/2) ■ ■ ■ upon making use of the periodicity of the 
integrands, hence the final result can be obtained. (Further details arc in Appendix C.) 

Moving on, we compute H ★ gk and * H . The first of these is easily seen to be 



H-kg n = n g v 



(72) 



but the other has an inhomogeneity. 



2nh n _i for n € N a 



dd 



g n *H = n g n + < 



(73) 



for n S N cvcn 

where the RHS involves non-diagonal versions of hybrid WFs, of the type indicated (see Appendix D and 
set 4>i {x) = expilx). Explicitly 



1 /" 27r 1 00 

hn,l (x, P) = Y ^(x-y) e~ a ^ e 2 ^dy = -e 2 ^~ l > ^ 
" j7r Jo k=0 



(-if f 
^k\{k + n)\ b2p - n+l+2k 



a 2i(p-l)x_ 



i-ir 



(74) 



2^p-i{p- i±a)!(p_i=a)] 

for integer (p — | (I + n)J > 0, but zero otherwise, upon using the series ©. We note that the right 
action of H has spread the g n distributions in momentum, through the effects of the inhomogeneities, and 
in particular has produced a contribution at p = — 1 from g$. 

In similar language, we could have abbreviated the explicit results in (|48|) and (|49|) as 



2n h-i. n for n£N 



del 



H*f n = n 2 f n + I 



(75) 



£n /l-2.n for n e 



2n „ for n S N 



dd 



f n *H = n 2 f n + < 



(76) 



£n /i-2.n for n e N c 
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where again the RHSs involve non-diagonal versions of hybrid WFs, of the type indicated (again see Appendix 
D and set <pi (x) = exp ilx) . Explicitly, upon using the series ([5]) , 

h^n (x,p) = ±- e^-y> Xn (x + y) e*»dy = 2 n ne 2i ^ x ^ ~*~ 1)1 S 2pil+n . 2k 

= 2 2 P -i ne 2Ki- P ) x \P2 — _U: (77) 

for integer (I + n) — pj > and (i (I + n) — p) < [n/2\ , but zero otherwise. 

It is significant that the RHSs of (|73[) . ([75)l . and (|75|) only have support for p <E Z, and so conform 
to our original expectation about momentum quantization. But in general, (|74p and (|77[) do not vanish 
when 2p £ Z, where contributions at semi-integer p can occur when I + n is an orfc? integer. So, upon 
considering these more general, nondiagonal WFs, the phase space must be expanded to include all points 
(x,p) G S 1 x Z/2. 

7.2 More hybrid WFs 

As anticipated in our abbreviated expressions for the various inhomogeneities, ()74j) and (|77|) . we may continue 
the process of hybridizing WFs by constructing Wigner transforms of other pairs of functions. If one function 
is an imaginary Liouvillc cigenfunction while the other is a free particle solution, we are led to define the 
hybrid 

1 r 2n 

h n (x,p) = — J i> n (x - y) <f> n (x + y) e 2lw dy (78) 

where xj) n are the Bessels, and </>„ are free solutions, 27r-periodic, and orthonormal in the usual sense. 

1 f 2 * 

- ^2^™ ( x ) = n <t>n ( x ) ) ^ J ^ k ^ ^" ^ dx = Sk ' n 
Then by standard techniques we have 

H-kh n (x,p) — n 2 h n (x, p) = h n (x, p) * p 2 (80) 



h n (x,p)*H = n h n (x,p)=p *h n (x,p) (81) 

as well as 



h k (x,p)* h n (x,p) = f n (x, p) S k ,n (82) 
and so on. However, since the Bessels on the circle are not orthonormal, we do not have similar simple 



results for h n (x,p) * hk (x,p), although this does not really seem to matter in practice. 
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These * product results can be established without explicit forms for either the h n or the /„, but it is 
perhaps useful to have such expressions in hand. Unlike the situation for the /„, as given in (|12[) . there is 
more freedom in the construction of the h n since we have a choice between taking <j) n (x) to be e mx or e~ mx , 
for n > 0, or some linear combination of the two. Explicit results depend on how we make this choice. 
Perhaps the simplest choice, conceptually, is to take only right- or left-moving plane waves for the <p n . For 
example, when <f> n is right-moving, then <p n (x) = e~ mx for n > is left-moving, and 

hi (x,p) = —J J n (e^) e -*»<*+i0 e 2i ™dy = | >_ ./ ^^ for p > n (83) 

but vanishes for p < n. On the other hand, if we take only left-moving plane waves for the 4> n , then 
4> n (x) = e mx for n > is right-moving, and 

1 /-27T f— IV 

/ T ( MX— v) \ „in(x+y) „2iypj„._ V x 7 



h^(x,p) = - J n (e«-">) e'"(-+»> e M "»dy = I J - e 2 ^+">* for p > (84) 



/ V / a 4P2>! (p + n) 

but vanishes for p < 0. 

In conjunction with the h n (x,p), we also define the duals 

/f.2-77 
o 

where %„ are the duals to the ip n . Various star products follow immediately. 



h n (x, p) = ±- I <t> n (x-y) X n (x + y) e 2lyp dy (85) 



1 r 27T e n e 2l ( x+y ^ for even n > 

h n (x, p)*H = n 2 h n (x,p) + — (j> n {x-y) e 2lyp { (86) 

K 7 ° 2ne i ^ x+y '> for odd n > 



p 2 * /i„ (x, p) = n 2 h n (x, p) (87) 



hk(x,p)* f n (x,p) = h n (x,p) S k ,n (88) 



fk(x,p)*h n (x,p) = h n (x,p) Sk, n (89) 



h k (x,p)*h n (x,p) = g n (x,p) 5 k ,n (90) 



h k (x,p)*h n (x,p) = f n (x,p) 4,n (91) 



h k (x 1 p)^h n (x,p) = e n (x,p) S k!n (92) 
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where the last of these brings us back to the plane wave WFs we introduced initially, in ©. Also, all the 
h n and h n are respectively left- and right-*-mapped by the dual metric into the h n and their conjugates. 



R(x,p)*h n (x,p) = e n h n {x,p) (93) 



h n (x,p)-kR(x,p) = e n h n {x,p) (94) 



This follows immediately from (|89|) . (|88|) . and 

More explicit results again depend on how we choose the free particle solutions. Choosing right-moving 
plane waves, with n > as above, we have 

K (x,p) ee 1- e^-ti Xn (x + y) e 2i ™dy = 2 ^~ n ne 2l ^ x ^% (95) 
2tt J q [n-p)\ 

for < n — p < [n/2\ , but zero otherwise. Choosing left-moving plane waves, with n > , we have 

h^(x,p) = — [ e - in( - x - y ) Xn (x + y) e 2iyp dy = 2 2 P+ n ne- 2lx( - n+ rt (— + P Z ^ (96) 



2* Jo "" v ° 

for < — p < [n/2\, but zero otherwise. For these specific choices we compute the right action of H to be 

£„e^" +2 ^^2j3 n-2 for even n > 
h%{x,p)*H = n 2 h$(x,p)+ { (97) 

2ne i{ - n+1 >8 2 p,n-i for odd n > 



e„e l( n+2)x 5 2p ,-n-2 for even n > 
/ i L( x , J5 )*F = „ 2 ^ (X;p) + < ; (98) 

2ne^~ n+1)x 82 V ,-„-i for odd n > 



The inhomogencities here can be identified with nondiagonal free particle WFs, e 
& x . It is significant that the inhomogeneities have support only for integer p 



i & and e L 2 

But note that this 



10 



includes, e.g., p = —1 when n — even though only /i^' L = 0) 7^ 0. Once again the star product of H 
with the h n (x,p) has spread the distributions on the phase space to give contributions outside their initial 
momentum support through the effects of the inhomogeneities. (For h„ >0 {x,p) cases, the initial support 
was only for negative momentum, and the star action of H also gives only negative momentum contributions, 
but spread out nonetheless.) 



10 Although in general the nondiagonal free WFs have support at semi-integer p when i+n is odd: l (x, p) = e*^ n— ^ x S2p,n+l 
and (x,p) = e i ^~ n ^' x S2 Pf — n —l- So once again, considering these more general, nondiagonal WFs, the phase space must 
be expanded to include all points (x,p) £ S 1 X Z/2. 
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8 Direct solutions of the dual metric equation 

Basic solutions to (|4"2")l . or cquivalcntly (j3"T]) . are obtained by separation of variables. We find two classes of 
solutions. The first of these is non-singular for all real p, although there arc zeroes for negative integer p. 

R(x,p;s) = sPT( [ +p j ex P ^-^scos2a;j (99) 

For real s this is real and positive definite on the positive momentum half-line. The other class of solutions 
has poles for all positive integer p. 

Rather (x,p; s) = — exp Qs cos 2x^j (100) 

For later use we also compute the left- and right-actions of H and H on R (x,p; s). 

H *R(x,p; s) = (p 2 - ^s 2 + (p - i ) s cos 2a; + ^s 2 cos4x ) R (x,p; s) = R (x,p; s) *~H (101) 



2, 

Linear combinations of (f9"9"| and/or (|100p are also solutions. This permits us to build a "composite" 

metric from members of the first class by using a contour integral representation. For t > 

1 / >(0+) ~ e tT 
R(x,p;s,t) = — R(x,p;st) — dr (102) 

2tti J_ 00 t 

The contour begins at -co, with argr = — tt, proceeds below the real r axis towards the origin, loops in the 
positive, counterclockwise sense around the origin (hence the (0+) notation), and then continues above the 
real r axis back to — oo, with argr = +7r. By construction, R (x,p; s, t) actually depends only on the ratio 
t/s. Evaluation of the contour integral gives 

t , . V i 



R(x,p;s,t)= --icos2x ^ (103) 

V s / (r(i+p)) 2 



where we have made use of 

1 1 



r(l+p) 2iri 

From R (x,p; s, t) we therefore recover our original dual metric by setting s = 2t. 



(o+) 

T- p - 1 e T dr (104) 



~ , v ( sm2 x) P ~,s , 
R(x,p:2t,t) = -i '—*=R(x,p) 105) 

9 The it root of the dual metric 

9.1 S as a direct solution of an entwining equation 

We look for an equivalence between the Liouvillc H = p 2 + e 2tx and the free particle H = p 2 as given by 
solutions of 

Six^y 1 *H*S{x,p) =p 2 (106) 
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or, barring complete invertibility, as solutions of the entwining equation (cf. (|159[) below) 

H * S(x,p) =S (x, p) * p 2 (107) 

This is again a special case of the ultra- local two star equation as given in [7] . For the Liouville-free-particle 
case, this amounts to an equation similar to that for R, but inherently complex. 

2ipd x S(x,p) =e 2lx e- 9 "S(x,p) = e 2rx S(x,p - 1) (108) 

Once again solutions are easily found through the use of a product ansatz. For any value of a parameter s, 
we find two immediate solutions: 

S (x,p;s) = sPT( [ +p ^ CX P (-\sexp(2ix)j (109) 

Mother (x,p; s) = -^r (— p) cxp Qscxp {2ix)j (110) 

The first of these is a "good" solution for p G (— l,oo), say, while the second is good for p 6 (— oo,0), 
thereby providing an overlapping pair of solutions that cover the entire real p axis, but not smoothly or 
continuously. These solutions could always be multiplied by periodic functions of p, exp (2iimp) for n G Z, 
but for integer- valued p (which is our primary interest) this has no effect. 

The two solutions for S are brought closer in appearance by shifting s — » sexp(±i7r/2) and using the 

— 7T 

reflection relation for the T: T (—p) = — ; r Thus we may take as our two basic solutions the 

r (1 +p) sin7rp 

alternate forms 

S (x,p; s exp (+*7r/2)) = e^ 2 ^ * - - exp (-is exp (2ix)J (111) 

-Sothcr (x, p;s exp (-nr/2)) = e^ 1+ ^ 2 (— — ) ) r- exp ( —sexp (2ix) ) (112) 
7r \siinrp J s p l (1+p) \ 4 J 

= e ^(i+2p)/2 A _M S(x,p;sexp{+iiT/2)) (113) 
\ sm np J 

which now coincide at the point p = —1/2, since we have also adjusted the normalization and phase of 
the 2nd solution. However, their derivatives with respect to p still do not match at p = —1/2, due to 

d_ (_ e ^(l+2 P )/2 1 sin ^ | ^ = i?r _ 

Taking all this into account, and exploiting the linearity of the equation, (|108p . a more general solution 
would be 

S—»i(*,p)=E J ds (c» 00 e-^ 4 + d n (s) e^ 4 ^) x e 2 "" ^ ^ exp (--^cxp (2fa)) 

(114) 

For integer-valued p the sum over n may be omitted. 
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9.2 The dual metric as an absolute it square 

Each such solution for S leads to a candidate real metric, given by 

R = S*S (115) 
To verify this, we note that the entwining equation for S, and its conjugate S, 

H*S = S(x lP )*p 2 , p 2 *S = S*H (116) 
may be combined with the associativity of the star product to obtain 

H*S*S = S{x,p)*p 2 *S = S*S*H (117) 

Thus the form in (|115p yields a solution to (|4"2")l . We need to work out a representative S * S star product 
to understand the relation to the previous solutions for R. 

We do this for the first form of the basic S solutions, (|109p . using the standard integral representation 
(|104[) for 1/r. We find a result that coincides with one of the composite dual metrics (|103[) . 

~ / 1 1 V 1 

S (x,p; s) ★5 (x,p; s) = I 1 H — -s 4 — —s 2 cos 2x 



16 2 J s 2p(T(l+p)) 2 

1 

1? 



R(x,p;s 2 ,l + -!-s 4 ) (118) 



By choosing s = ±2, we again obtain the original dual metric. 

S (x,p; ±2) * S (x,p; ±2) = R (x,p; 4, 2) = -A- > = R (x,p) (119) 

(r (p + 1)) 

This provides a greater appreciation of the information contained in S, and motivates us to consider an 
alternative construction of such roots of R. 



9.3 S as a sum of hybrid WFs 

The basic ideas here are essentially the same as used in the construction of R as a sum of f n , only the 
bilinears appearing in the Wigner transforms involve two different types of functions: One wave function is 
an imaginary Liouville cigcnfunction while the other is a free particle solution, precisely the h n (x,p) defined 
earlier in ([78]) . 

We form the sums 



S(x,p) = ^2 sfe^i h n (x,p) , S(x,p) = ^2 h n (x,p) (120) 



n=0 n=0 
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Then from ||5D|) and l|8ip we immediately obtain 



H-kS(x,p) = S(x,p)-kp 2 , p 2 -kS(x,p) = S(x,p)-kH (121) 
and, upon using 11820 . we also obtain R as the absolute star-square of S. 



S (x,p)*S(x,p) = ^2 V^™ h n (x,p)*^2V~£k h k (x,p) = ^2 £ n fn {x,p) = R (x,p) (122) 



n=0 k=0 n=0 



However, since the Bcssels on the circle are not orthonormal, we do not have a similar relation for S (x,p) * 
S (x,p). While these results can be established without explicit forms for the h n (x,p) and for S {x,p), it is 
perhaps useful to have such expressions in hand. 

Like the individual hybrid WFs, explicit results for the sum depend on how we choose the free particle 
solutions. In particular, we may take the hybrid WFs built from right-moving and left-moving plane waves 
in (|53")) and (|84[) . This gives rise to the phase-space equivalent of the "chiral kernel" in [8]. 

S R (x,p) = £ K (x,p) = UL e 2 ^ yfr - } ■ ^ 2 " lX 



4?p\ ^—i (p-n)\ 

n=0 1 n=0 yl ' 



= (l-V2) ^e 2,pl + yfi—^—te-^T (p+ l,--e 2ix ) (123) 
where the incomplete T function has made an appearance (e.g. see [1]). 

poo z p+l 

T{p+l,z)= S p e- s d S = T(p+l) lFl (p+l;p + 2;-z) (125) 

We leave it as an exercise for the reader to use the properties of the incomplete T to check directly that 



Sr{x,p)* S R (x,p) = R(x,p). 
Similarly, we have 



~ (—1) 1 

s L (x, P ) = £ ft£ = [ -^zr^ px E ^ 2 ^h + nv e2mx (126) 

n=Q y ' n=0 KP >' 

-tilled + ^izilLe^" ( Tip+1) -r(p+l,le 2 ™)) (127) 



2 



4p (p!f 2p (p\) 

Comparison with the results in [8], §4, Eqn's (67) and (69), shows that Snj J (x,p) are essentially Fourier 
transforms with respect to one variable of a particular combination of Lommel's functions of two variables. 
Who knew? 

Alternatively, we could take the free particle solutions to be non-chiral: 1 for n = but for n > 0, 
(e-" lx + (-1)™ e mx ) /\/2. That is 

K (x) = ^ (e~ mx + (-1)" e mx ) (128) 
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which are properly normalized, as in (|79p . In this case we would be led to the Wigner transform of the 
well-known generating function. 



oo 

Az sin x 



Jn(z)e mx = J a (z) + J2 Jn(z)(e mx + (-l) n e- mx ) (129) 

71 — — OO 71 — 1 

OO 

= V^Jn (Z) K (X) (130) 

n=0 

Thus we obtain a remarkably simple result for the non-chiral kernel. 

S NC (x,p) ee -L f 2 J e^ v) sM*+y) e 2*Py dy = IxL exp (131) 
which wc recognize as a particular case of our basic solutions, (|109l) . 

S NC (x,p)=S{x,p;s = -2) (132) 
The selected non-chiral free particle solutions have WFs with symmetric momentum support. 

e„ (x,p) = (5 p>n + 5 p _ n + 2 (-1)" cos (2nx) x S pfi ) (133) 



The corresponding hybrids arc linear combinations of those in (|83|) and (|84|) , as given by 
C (x,p) = ^ K (x,p) + (-!)" ft* ■(*,„)) 

_I^_^ e 2 l ( P+n ), for n>p >Q 



(134) 



(~ 1 ) P e 2i(p+n) a; , 2"(~l) p " + e U{p-n)x r > 

^ 4f2™p!(p+n)! e + V ^ 4Pp!(p-n)! +e 101 - " 

and vanish for p < 0. All this leads us once again to the expected relation (|122[) . 



10 Meanwhile, back at the metric 

We might also wish to systematically determine all metrics, rather than dual metrics, through which the 
Liouville Hamiltonian is rendered hermitian. This is perhaps a more conventional problem to attack in the 
framework of quasi- hermitian theories [TBI [IH E] . That is to say, we seek all solutions to 

H*R = R*H (135) 



for real functions R{x,p) = R(x,p). (This should be compared to (14211 given above.) In this context it 
is somewhat repetitious but perhaps instructive to go through a few details omitted in the previous two 
sections, including some forays into calculational cul-de-sacs. 
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10.1 R as a formal sum of the dual WFs (Not!) 

In parallel to the previous construction of the dual metric as a sum of WFs, (f3"4")l . we might try 

oo 

Q(x,P) = ^£klk (x,p) (136) 
which gives, at least formally, the expected phase-space orthogonality relation 



rt,Q(x,p)fn(x, P ) (137) 



However, Q does not satisfy the homogeneous equation (|135|) due to the inhomogeneitics resulting when H 
and H act on the individual f^. Moreover, the best we can do with the sum (|136p . so far, is to interpret it as 
an asymptotic series related to the Wigner transform of an integral whose asymptotic expansion is a formal 
generating function for the dual wave functions. Rather than pursue this here, we turn to direct solutions 
of the differential-difference equation corresponding to (| 135[) . 

10.2 Solving directly for R 

We find the following left- and right-sided star products, and their difference. 

2 | -2ix\ , D / \ _ I I „ * q \ -2ix„d„ 



H*R(x,p) = ( P ' + e-* x )*R{x,p)=Up--d x ) +e- Mx e°*jR{x,p) (138) 

R(x,p)*H - H*R(x,p) = 2ipd x R (x,p) + 2i sin (2x) R(x,p+1) (140) 

So then, solving (|135[) amounts to solving the linear differential-difference equation^ 

- pd x R {x,p)= sin (2x)R(x,p+l) (141) 

This should be compared to (|3T|) given above, which it becomes upon letting p — > —p. (Note that R actually 
corresponds to R^ 1 in that earlier discussion.) This is immediately solved upon assuming a product form, 
R (x,p) = q {x) r (p). There are two distinct sets of solutions corresponding to positive and negative constants 
of separation, ±s. 

The first one-parameter (s) set of solutions is 

R(x,p;s) = s p T(p)cxp (^scos(2x)J (142) 



11 This first-order equation becomes a second-order equation if the Voros product is used. That is, if we demand H<\R = R<H 
then we have to solve: (9^ — 2ipd x — exp (2ix)) R (x,p) + exp (— lix) R (x,p + 2) = 0. Moreover, the solutions of this second- 
order equation are not real. So we prefer to use * and not <. 
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up to an overall multiplicative constant. Identification with the previous dual metric basic solution is given 
by -Rothcr (x, — p; s) = R {x,p\ s). Also, the duplication formula V (p) = (^p) r (| + ^p) gives 

R (x,p\ s) = ^jLr (±p) r (i + \p) exp Qs cos (2a:)) (1 43) 

for whatever that's worth 12 while the reflection relation T (p) — ^ — 2_ gives 

1 — i \f / "(1— p) sin Trp ° 

7TS P /l \ 

i? (x,p; s) = — — exp -s cos (2a;) (144) 

1 (1 — p) smnp \2 ) 

For real s, the solutions (|142[) are hcrmitian and positive-definite functions of real variables x,p so long as 
p > 0. However, R(x,p;s) is not bounded either on the negative p half-line, or on the positive p half-line 
for any real s (cf. Stirling's approximation). Otherwise, (— l) fe+1 R(x,p) is hermitian and positive-definite 
for — k — 1 < p < —k, where k = 0, 1,2, • • ■ . 

Now we know from previous considerations that for n > the periodic WFs f n have support only for 
positive p, so (|142p would seem to be a preferred set of solutions, except for the ground state: fo has support 
at p = 0. This issue must still be addressed. 

Another set of solutions is given by the form 

s p ( 1 \ 

Mother (x,p;s) = —r exp ( --scos(2x) 1 (145) 

Identification with the previous dual metric basic solution is given by R(x, —p; s) = i? thcr (x,p;s). The 

duplication — ; = —r, — , 2 w — i — 7 and reflection relations — ; = r (p) sln7r P now m VC 

y T(l-p) r(i-ip)r(i-i P ) r(l-p) T 8 

Mother {x,p; s) = —jj — ^\ S \ yf — exp (^-iscos (2x)] (146) 
_ , . s p r (p) sin 7rp / 1 \ 

Mother (a;, p; s) = exp I - - s cos (2x) I (147) 

again for whatever that's worth. For real s, these other solutions are hermitian and positive-definite functions 
of real variables x,p so long as p < 1. Otherwise, (— l) k i? thcr {x,p; s) is hermitian and positive-definite for 
k < p < k + 1, where k = 1, 2, • ■ ■ . Also, i? ther is now bounded for all real p. 

An interesting problem now is to find the real square-root of either R (x,p; s), or i? thcr (x,p; s). This is 
not just s p / 2 y^T (jp)eT scos ( 2x \ say, since the square-root must be taken in a * sense. That is, we seek a real 
S such that 

R(x,p;s) = S (x,p;s)* S (x,p;s) (148) 



12 These solutions immediately call to mind the Mellin transforms of Besscl functions, J v (p) = r f — r , which arc 

' y " 2 pr(i+^) 

solutions to the Liouville energy eigenvalue problem in momentum space, (p 2 + e~ 2S p) J v (p) = v 2 J v (p). 
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Were it not for the T function, this would be easy since 

exp (a cos (2») + bp) = exp ( 2co ^ (fc/2) cos (Jte) + * cxp ( 2co ^ (fc/2) cos (2*) + (149) 

But the presence of the T function makes the problem a little more challenging. 

Recall for p > 0, T (p) = J °° t p ~ 1 e~ t dt while for all p there is again the contour integral representation 

cmi) , so 

1 /-(o+) 

i p_ Vdt (150) 



r(l-p) 2?ri 

The contour here is the same as in (|102|) . Thus 

/■(0+) ( l \ 

Mother (x, p; s) = J t p 1 cxp It + ps - -e s cos (2x) J dt 

r(o+) / i \ 

= j cxp ( e) nt + pint + ps - -e s cos (2a;) J dint (151) 

Star composition of two such integrands is now possible, as in (|149[) . although we then have to deal with a 
double integral.. But taking the square-root in this approach is not transparent. 
From another perspective, the problem involves computation of either 



/>oo />oo 

r(p)*/(z)=/ t p ~ 1 -kf(x)e t dt= j exp (pint) */ (x) e* dint 
Jo Jo 



J exp (pint) f ^x - ln/j e*d Int = J 



t p ~ 1 e t f[x-h,lnt 



-*f(x)= t p - 1 *f(x)e t dt= exp (pint) * f(x)e t d Int 

r(i-p) J-oo J-oo 



r(0+) / i \ r(0+) 

J cxp (pint) / ( x — — i Int J e dint = J 



1 



t p - 1 e t f[x--i\nt) dt (153) 



Alternatively, since V (p) * e lkx = e lkx T (p + i/c), then if / (x) = J e lkx F (k) dk the two computations become 
r (p) * f (x) = j Y (p) * e lkx F (k) dk = J e lkx T (p + ^kj F (k) dk (154) 

*f(x)= f — — ^ * e lkx F (k) dk = f e lkx - F( - k \ dk (155) 



r(i-p) JKJ J r(i-p) ' w J' r(i-p-|fe) 

We also note that 



exp (a cos 2x) = V I n (a) e 2tnx , I n (a) = - / e a 

K /n 



cos 



cos (n6) d0 (156) 



n— — oo 



with (a) = /_„ (a) and 7„ (—a) = (—1)™ I n (&)■ Therefore we may write, as least formally, the ordinary 
product in R as a star product 

oo oo -p / \ oo T ( \ 

r (p) exp (a cos 2x) = £ /„ (a) T (p) e 2m * = £ 7 n (a) - ^ T (p) * e 2m * = T (p) * ]T e 2i " x 

71 — — OO 71 — — OG ^ ' 71 — — OO 

(157) 
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where the Pochhammer symbol is (p) n = (p) (p + 1) ■ • ■ (p + n — 1) = T (p + n) /T (p). Note the zeroes in 
the summand for integer p > when n < —p. Unfortunately, this is not of much use without a closed form 
for the sum. This would appear to be one of the aforementioned cul-de-sacs, so let us try a different route, 
parallel to that used previously to construct S. 

10.3 Solving for VM 

We look for an equivalence between the Liouville H = p 2 + e 2lx and the free particle H = p 2 as given by 

H = S+H+S- 1 (158) 

That is to say, in phase-space we wish to solve the linear equation (cf. (|107|) above) 

p 2 * S = S * H (159) 

This is again a special case of the ultra- local two star equation as given in [7] . For the Liouville-free-particle 
case, this again amounts to an equation similar to that for R. but inherently complex. 

-2ipd x S(x,p) = e 2lx e d ?S(x,p) = e 2lx S (x,p + 1) (160) 

which becomes the previous (|108[) upon p — > —p. Once again solutions are easily found through the use of 
a product ansatz. For any value of the parameter s we have two immediate solutions: 

S 0, p; s) = s p T (p) cxp (-s exp (2ix) ) (161) 



Sothor {x,p; s) = -yexp (^--sexp (2ix)j (162) 

The first of these is a good solution for p 6 (0, oo), say, while the second is good for p 6 (— oo, 1), thereby 
providing a solution for the entire real p axis, albeit a discontinuous one. The relation to the previous dual 
solutions is obviously 

S(x,p; s) = 5 thcr (x, -p; s) , S otilCY (x,p;s) = S(x,-p;s) (163) 

As in the dual situation, these solutions could always be multiplied by periodic functions of p, exp (2iTmp) 
for neZ, but again, for integer- valued p this has no effect. 

The two solutions for S are brought closer in appearance by shifting s — > se ±l7T / 2 and using the reflection 
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relation for the T. Thus we may take as our two basic solutions the alternate forms 



S{x,p;s exp (+wr/2)) = e l7Tp/2 s p T (p) exp f^s cxp (2ix)j (164) 



pother (a;,p;sexp(-wr/2)) = e ra(1 " p)/2 s p sin (tip) T (p) cxp Qsexp(W)j (165) 

= e*^ 1 -*)/ 3 sin (Tip) 5 (x,p; sexp (+»tt/2)) (166) 

which now coincide at the point p = 1/2, since we have also rescaled the 2nd solution. However, their 
derivatives with respect to p still do not match at p = 1/2, due to ^ (e™' 1-2 ^/ 2 sin (717?)) 
Taking all this into account, and exploiting the linearity of the equation, a more general solution would be 

general (x, P ) = ]T J d S (c n (s) e"^ + d n ( S ) e^ 1 "^ 2 sin (^p)) X e 2 ™* S*T (p) 6X P ^SCXp (2ix)j 

(167) 

Once again, for integer-valued p the sum over n may be omitted. 

Since these solutions for S are complex, a suitable real metric would now be given by 



= —ITT. 

p=l/2 



R = S*S (168) 

and as before we have H*R = R*H consistent with S-kH-kS^ 1 = H = H = S~ 1 -kH*S. For completeness, 
we extend this to a simple theorem that can be used to construct other composite solutions for the metric 
through multiple star products: 

[Lemma] Any solution of (|159p , i.e. p60[) , gives a real solution of (|135[) . i.e. (|141[) . namely R = S*S. 
Moreover, appropriately ordered odd star products of any solutions to (|135|) . (|159[) . and the complementary 
equation (|42|) are also solutions to (| 141 [) and (| 1 60[) . For example, H * i?i ★ R2 * R3 = R\ * R2 * R 3 * H. 
Similarly, H * Si * R 2 * R3 = Si * R 2 * R3 * H . 

We need to work out a star product at least for a particular case just to understand the relation to the 
previous solutions for R. So we do this for the first form of the solutions, (|16ip . using the standard integral 
representation for T (p) . We find a composite metric 

R = S (x,p; s) ★ S (x,p; s) 

i_I s 2 e 2 fa ) p ^ _ I s 2 e - 2 ^ P s 2p r 2 Cp) 

= (l + ^s 4 - is 2 cos 2x\ " s 2p T 2 (p) (169) 

assuming s is real and provided Re (s 2 e ±2lx ) < 4 to avoid the singularities in the (l — j.s 2 e ±2tx ) P factors. 
For example, if e 2s < 4, then all real x satisfy this condition, and the composite metric is real and positive 
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definite for real p > 0. The double pole at p = is somewhat troubling, and again presents an issue for the 
ground state WF. Nevertheless, we easily check that this composite R satisfies the metric equation (|141|) . 
Similarly, for the other solution (|162p . we find 

R = Mother {X,p) * SUhcr (x,p) 

1 \ ~ p / 1 \ _p s 2 P 



4" " J y 4~ ~ J T 2 {l-p) 

1+ ir-2 sc H jvjtt) (170) 

again for real s and provided Re (e 2s e ) < 4. Again we easily verify that this composite R satisfies the 
metric equation (|141| . We also note that these two composite metrics are more simply related than the 
previous forms of the metric or even the pairs of solutions for S. Namely, we go from one composite R to 
the other just by the interchange T 2 (p) <-> 1/r 2 (1 —p). 

We should also express the composite solution in terms of (sums) of the previous solutions for R, if that 
is possible. Indeed, just as in the dual metric situation, it is possible. We find 

^l + -L s 4 -i s 2 cos2x) P s 2p T 2 (p) = J™ ^e-( 1+ ^) T R(x,p; S 2 T) (171) 

where R (x,p; s 2 t) is as defined in (|142[) . Also 

/I 1 \ ~ p s 2 P / >(0+) dr i i 4i 

( 1 + i6 2 C0S2X ) T 2 (l-p) = / x i^ +iss)T Ro^{x,p-s 2 r) (172) 
where i? ther (x,p; s 2 t) is as defined in (|145|) . 

10.4 VR as a sum of hybrid WFs (Not!) 

In parallel to the previous construction of S as a sum of the hybrid WFs h n , as in (|120[) . we might try to 
construct roots of R as formal sums of the dual hybrids h n . 



T (x,p) = ^2 Vs7i h n (x,p) , T(x,p) = ^2 h n (x,p) (173) 

n=0 n=0 

However, these are not solutions to (|159[) or its conjugate due to the inhomogeneities resulting when H and 
H act on the individual h n and h n . Even so, these and similar sums do yield some interesting star products, 
at least formally. For example, 

oo 

T (x, p) * T (x, p) = J2 eJk (x, p)=Q(x,p) (174) 

k=0 

Rather than pursue this here, we consider other direct solutions of the differential-difference equation corre- 
sponding to (|159D . 
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10.5 Additional solutions for S 



Rc-instating the coupling constant m via exp (ix) — > mexp (ix), as in ([T]), a straightforward series solution 
in powers of m gives another form for solutions to (|159|) . 



S(m;x,p)=T(p)(J-^ e"^% {me") 



(175) 



S othcr (m;x,p) = — - - (-) (-if e -*% (ime lx ) (176) 

1 (1 — p) \?71 J V ' 

the latter obtained from the first form by x — » a; + 7r / '2 and r (p) — > 1/r (1 — p). We also note that (|176D 
may be written as 



5*oth or (m;x,p) 



1 



2\ p 



"J v (me 1 



r (1 — p) \m / 

Here we have made use of the modified Bessel function l v (z) with properties 



d 7 ( 

z—I v (z 

az 



In {z 



exp (zcosf 



(. 2 /4) fc 

\2J ?-jkW(v + k+l) 



k=0 



(-z 2 /iY 



= { -iyi u (iz) = (-) £ fc!I> + fc + 1) 



= zi^ii (z) ± ^ (z) 



1 



2 71 



= — / e zcose e ±iTie d6» for n £ N 

oo 

= /„(z) + 2^/ t (z)cos(M) 
fe=i 



(177) 

(178a) 

(178b) 

(178c) 
(178d) 
(178e) 



Either the first or the third of these properties leads to the most direct verification that (|175| and (|176l) 
satisfy the equation (|159[) . 

We also recall the elegant contour integral representations of Schlafli and Sonine which are valid for all 
v and z. 



Iu{z) 
Ju(z) 



zy l f {0+) _ v _ x ( 
— - — / w exp w + - — aw 



z\ v 1 



(0+) 



w " " exp [ w — - — I dw 

4w , 



(179) 
(180) 



.2/ 2m 

The contour is the same as that used in (|102|) . Thus we identify (|175| and (|177j) as linear combinations of 
the basic solutions (|161[) and (|162[) (after re-instating the coupling constant m via x — > x — i lnm), expressed 
as integrals over w = l/s. 



S (m; x,p) 



1 

2?ri 



1 r(o+) 
Mother (m;x,p) = — / S, 

27T* J-oo 



{0+) q f n > 

o x — « m to, p; s = — e — 
-oo V wj w 

other x - zlnm,p;s = — e 



(181) 
(182) 
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In the same manner, we can construct additional solutions to (|107|) , with m re- instated, to obtain S (m; x, p) = 
Mother (m;x, —p) and Sother (m;x,p) = S (m;x, —p). We leave as an exercise the calculation of the corre- 



sponding composite metric, R(m;x,p) = S (m;x,p) * S (m;x,p). 

11 Metric dependence of expectation values 

It is important to realize that expectation values are metric-dependent, in general. As a first illustration of 
this, we consider expectations of WFs for energy eigenstates using the basic solutions (|142p for the phase- 
space metric. 

M„ 0) = 2^:^ R(x,p;s)f n (x,p) = -^^2 sPT (p)J fn(x,p)exp 0scos2xJ dx (183) 

For the ground state case, n = 0, there is a problem here. The sum over p must include p = and this 
diverges due to the T (j>) in R(x,p; s). The individual basic metric solutions are not useful in this one 
case. It would appear to be necessary to take a linear combination of the basic metrics an integral over 
s, perhaps - to tame this singularity. We leave this as an unsolved problem. 

For all n > 0, the individual R (x,p; s) give finite, reasonably well-behaved A^n (s). We have 

Mi (s) = it R(x,p;s)f n (x,p) = jTs p T(p) j f n (x,p)exp C^scos2x) dx (184) 

From the series (|12l) and the integral representation of the modified Bcssels, (|178d[) . these basic norms are 



°° s p F (p) p ~ n 1 /s\ 

Nn (S) = S ~^r- E k \ { n + k)\(p-k)\(p-k-n)\ X In - p+2k (2) 



(185) 



fc=0 

Alternatively, 

z sin 2a; \ 



^ s p 1 /isin2a;\ p /l \ 

M n (s) = > —7 rr — / ( ) LegcndrcP (p, — n, i cot 2x) exp I -s cos 2x ) c?x (186) 

^ n p(p-n)\2TT J \ 2 7 V 2 / 

1 °° 1 /1 \ 1 r2ir . \ / 1 \ 

= -7 777 tt I 7 s I 77- / e 2lfca; KummerM ( n + fc, 1 + n, ^e^ 2 " ) exp ( -s cos 2a? | dx 

nl kl (n + k) \A J 2tt J Q V '4 / V 2 / 

The infinite sums here are convergent, but we have not found a closed-form. Nevertheless, the norms are 
clearly s- and hence metric-dependent. If we switch our attention to the dual WFs, we encounter only finite 
sums of Bcsscl functions. 
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11.1 Diagonal dual WFs 

As further illustration of the metric dependence of phase-space expectations, we compute expectations of 
dual WFs for energy eigenstates using the basic solutions (|M|) for the phase-space dual metric. 

A^(s) = 7T~^ R{x,p;s) f n {x,p) = -^it, s py (i + p ) f fn (x, p) exp cos2x j dx (187) 
From the series (jTTJ) and the integral representation ()178d[) these dual norms are 

p=0 y 11 fe=max(0,n-p-Ln/2j) 

Note that A/" n (s) = M n (— s). For example: 



A/"o («) = /o (2 J (189) 
M (*) = -4 d) (190) 



a u V2, 

^ ( . ) _« ( ^ +8)Jb (i)_H Jl (') 



;|) + _L (8 ^ 2)/ „(£ 

12 /„ 9 > T /s\ 288 /s 



8^(2J+^( 8 - s2 ) J o(o) 



= £(3.» + 32) Jo(f)-^/i(! 



^Mi) + S( 32 - 3s2 ) /o (0 < 1M > 
These results should be compared to the much simpler norms (|33[) obtained upon computing phase-space 
averages using (|30[) . The basic norms are indeed positive for all s, but s dependent, as is evident upon 
graphing the chosen examples. 

11.2 Non-diagonal dual WFs 

We also compute, for k 7^ n, 



1 



A4, n (s)= — T R(x,p;s)f k . n (x,p) =Af n ,k(s) (193) 

27r Jx,p 

to see if they vanish. As a consequence of the inhomogeneities in the dual energy eigenvalue equations, it 
is not obvious that they should vanish for k ^= n. In fact, they do not. 

To verify this statement and to compute A4,n (s), for general k 7^ n, we first determine the non-diagonal 
dual WFs. 

fk,n (x,p) = — J Xk(x-y) Xn (x + y) e 2vyp dy 

Lfe/2J L«/2J . u . , , 

= tfW<*-0 J] J] (fe-J-l)-("-^-l)- e ^(i-J) Sk+n _ 2pt2j+2l (194) 
i=o z=o ■ 7 ' ' 
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The sums result from using the series in ([5]). Given the range of the sums, we must have < 2p < k + n 
for a non-zero result. So, as for the k = n case, the support in p is finite. Also, for a contribution 
to the sums, we must have k + n = 2pmod2. Thus fk. n (x,p) can be non-zero for semi-integer p, when 
k + n is an odd integer. Let us simplify the double sum in (|194|l . First, suppose k + n is even. Then 
Sk+ n -2p,2j+2i under the double sum implies p must be integer and < ^y 1 - The upper limits of the sums 
further imply k/2 — [k/2\ + n/2 — [n/2\ < p. That is to say, for k + n even, we must have p integer and 
< (k + n) /2 — \ k/2\ — \n/2\ < p < ^p 1 - On the other hand, for k + n odd, we must have p semi-integer 
and \ <p < -^Tp- . Incorporating these conditions and using the Kronecker delta to eliminate one sum, for 
both even and odd k + n cases, we find 

*- (*»&- W -L,j) f (» + L*5-J-w-J)i 

(195) 

Note that this coincides with (fl7)) when k = n. Also note we do not need the floor function in the [^p-J 
factors when k + n is even, but we do need it when k + n is odd. Similarly, we do not need the floor function 
in \p\ when k + n is even, but we do need it when k + n is odd. In particular 

f^i(x,p) = 2e- ix 5 P!l/2 (196a) 

fa{x,p) = 2<5 p ,o + 8e~ 2ix 5 pA (196b) 

fiZ (x,p) = 4e ix 5 pA/2 + l6e- lx 6 p , 3/2 (196c) 

fiZ 3 (x,p) = 126 pA + 96. 2 %, 2 (196d) 

^ = 12e-^5 Pil/2 + (48e" + 96 e - 3 ") 5 P)3 / 2 + 384e~ lx 6 p , 5/2 (196e) 

(a;,p) = 4<5 P , + (l6e 2 " + 64e" 2M ) S pA + (256 + 768e" 4M ) S p , 2 + 3072 e - 2M <5 p , 3 (196f) 

We now finish checking nondiagonal orthogonality for the basic metrics, i.e. whether A4,« (s) = for k =/= n. 
This is not true, in general. 

First, suppose k + n is even, and consider the simplest case. The average over position is 



fo, 2 {x,p)R (x,p; s) dx = — - — / (2S p . + 8e 2 "<5 p ,i) exp (- - cos 2x ) dx 
2tt J s p T (1 +p) 2n J \ 2 ) 

(2Io (|) S p . - 8/1 (|) (197) 



spr(i+p) 

This is real for real s, so 

1>27T -1 /"2-7T 



75-/ /o,2 {x,p) R{x,p;s)dx = — / ., 
2tt J 27r J 



f 2 ,o(x,p)R{x,p;s)dx (198) 
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Summing over p now gives 



M), 2 (s) = Af 2 , (s) = 2/ (I) -h = 2/ 2 f i ) (199) 



-2/ s V2/ V2 
This does not vanish for s ^ 0. Similarly for other even k + n cases: 

J=max(0 ,^_ p _ L fJ) 1 2 j 

We leave odd /c + n as an exercise for the interested reader. 

In contrast to (|200p . all the non-diagonal dual WFs vanish when summed/integrated over phase space 
using the dual metric of ([30)) . 

^-iR(x,p)f^n(x,p) = for k^n (201) 

as follows most easily from ([34]) . For the case when k + n is even, this involves cancellations among the 
various terms for different p after integrating over x. For example, 

I /-27T _ x |.2tt (sin 2 x) P 

o~ / y]/o. 2 (x,p)R(x,p)dx = — / / ,2 (a;,p) 2 cfa 

1 /" 27r 1 /" 27r 

= — / 2dx+ — Se~ 2 ' lx sin 2 xefe = 2 - 2 = (202) 
2* Jo 27r7 V 7 

However, when k + n is odd, and therefore p is semi-integer, we interpret the metric as integer powers of 

|sinx|. Then the individual terms contributed by fk t n (x,p) to the sum over p each vanish separately upon 

integration over x. In this regard we note that Tw^kp is still a solution of (|3ip . for p > 0, given that 

(s'mx) 2p 6 (x) = 0. (That is, all test functions are required to be non-singular at x = 0.) 

Thus it is evident that requiring the phase-space averages of fk t n to vanish for k ^= n imposes additional 

conditions on the dual metric, beyond those specified in (|31[) . However, these additional conditions are not 

sufficient to force the dual metric to be proportional to R(x,p) in (|30|) . Any dual metric of the form given 

in (|34p but with arbitrary coefficients, i.e. R(x,p) = X^fcLo c kfk (x,p), will give a vanishing phase-space 

average for nondiagonal dual WFs. This follows immediately from 

7T- p fj (x,p) fk,i (x,p) = S jjk Sj,l (203) 

From (|200|1 we can only conclude that R(x,p;s) cannot be so expressed as a linear combination of the 
diagonal Of course, this can also be established by other means. 

As a consequence of A4.n (s) ^ the norms of general pure state dual WFs are even more complicated. 
For example, consider a linear combination of dual functions, \ — aJ ^o + PA 2 , and form the corresponding 
dual WF. 

7=^-1 X(x -y) X{x + y) e^Pdy = \a\ 2 % + + ctffo + |/3| 2 h (204) 



2tt 



o 



3G 



Under the action of R (x,p; s) the normalization of this pure state dual WF is 

1 r 2 * ~ ~ ~ ~ _~ 

Af{s) = — R(x lP ; S )f(x,p)dx=\a\ 2 Mo(s) + af3No,2(s)+aPN2 fi {s) + \P\ 2 Af 2 (s) 
Zn 



u 



«r + ? (8 - s ) W) h (2 J + + 2«/3 + 8 |/3| Z J h J ( 205 ) 

On the other hand, under R(x,p) in (|5D|) . the norm of this dual WF would be just \a\ 2 + 2 |/3| 2 . 

Ultimately, A4,n 7^ for k ^= n originates in the inhomogencities in the dual eigenvalue equations. By 
contrast, a similar calculation for the metric, in those cases where is it is well-defined (k =/= 7^ n). would 
give 

A4, n 0) = ^- £ R (x, p; s) f k>n (x, p) = (206) 

for k ^= 11. This is a simple consequence of the Lone Star Lemma, the homogeneous equations H*fk t n (x, p) = 
k 2 fk,n {x,p) and fk,n (x,p) * H = n 2 fk,n (x,p), and (|135|) . Hence, for a linear combination of WFs corre- 
sponding to the wave function ip = aJ k + (3J n , i-e. 

f= h jj * {X -y)^'^ + y) e2lVPd V = M* A + aPfn,k + aM, n + l/f In (207) 
the action of R (x,p; s) would produce the normalization 

1 f 27r ~ 

AA(s) = — / R(x,p;s)f(x,p)dx= \a\ 2 M k (s) + \p\ 2 M n {s) (208) 





We have assumed that neither k nor n are zero. 



11.3 Expectations of H 

The corresponding expectations of H under the action of R (x,p; s) are given by 
tk,n {s) = 77- Y, R (x,p; s) (fk. n (x,p)* h) 

Z7T TZ „ V / 



K 

I x,p 

( 



— T, R(x,p;s) 
2lT Tx, P 



2n h-i,k for n E N dd 
n 2 /fc,« + { (209) 
£ n h-2,k for n E N 



where diagonal cases will be denoted by 7i nj „ = TL n . Alternatively, we may use the Lone Star Lemma and 
(fTUT]) to write 



Wfc,n (s) = ^ (h * R (x, p; s)j / fc ,„ (x, p) 



2tt ^p( P 8 



s 2 + ( p - i j s cos 2.t + is 2 cos4a; ) R(x,p;s), 



Since R{x,p]s) is real for real s, and = (|42]) then gives 



H fe ,„ (s) = W„, fc (s) (211) 
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For example, using (|196b[) we have 

~ 1 [ 27T ( 1 1 1 \ ~ 

^Uiiis) — — I --s 2 - -s cos 2x4- -s 2 cos4x )R(x,0;s) x 2dx 
27r Jo V 8 2 8 / 

1 f 27r /ll 1 \ ~ 

H / 1 - -s 2 + -scos2x+ -s 2 cos4a; ) R(x, l;s) x Se~ 2lx dx (212) 

2W V 8 2 8 J V ' ' ; V ; 

From (fl78d|> and (TiTSc]) we then obtain 

«... w = + r 2 Mf)) * 2 

i _ ' >) h (') + u m + u a) - -L.»i, fi) - 4a 3 m) x s 



2/ 4 U V27 4 V27 16 \2J 16 V2/y s 
8/2 (i) (213) 



So, for real s, 7io,2 (s) = Ho, 2 (s) = 7^2,0 (■?)■ Similarly, the n = and 2 diagonal cases are given by 

H (s) = ±s 2 I 2 Q (214) 
H 2 ( s )=48/ 2 (|)+^(8- S 2 )/o(|) (215) 

to be compared to J\f n (s) given above. The important point here is that the average of H for definite energy 
dual WFs, when computed as the ratio H. n (s) /Af n (s) using the basic metric R (x,p; s), is far from a simple 
expression. It is certainly not n 2 , as it was in (|52p . when computed using the dual metric of (|30l) . 
With these cases in hand, we can also consider the linear combination of dual WFs in (|204| . 

H{s) = ^^R(x,p;s) (7(a>,p)*Jj) = \a\ 2 H (s) + a/3H ,2 (*) + a/3W 2 ,o (s) + |/3| 2 H 2 (a) (216) 
We find 

H (s) = § (8 - S 2 ) |/3| 2 Jo + f ^ 2 |a| 2 + 8 (a/3 + a£) + 48 \{3\ 2 ) I 2 (1) (217) 



,2/ 

to be compared to Af (s) given above. Again, the important point is that the ratio Ti. (s) /TV (s) is far from 
a simple expression, and is not 4x2 |/3| 2 / (\a\ 2 + 2 \(3\ 2 ^j , as it would have been had we used the dual metric 
of ([llll- 
Other cases of even k + n are handled similarly, using (|195[) . Wc have 

J (p 2 - ^s 2 + (p-^j scos2x+ is 2 cos4.T^ R(x,p;s) f k , n (x,p)dx 

min(| % |.n+| I - LpJ ) , , x, 

2^ 



spr(i + ») ^ 7! (£±2 -p- j\\ 



X 2~ J if ~ + 2~) scos2x +l s2cosAx ^ e 2i{k ~ p ~ 23)x cxp ^-iscos2x^ (218) 
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Then (|178d[) and p78c[) permit evaluation of the spatial integral to obtain 



2jt 



271 



P 2 - + I P 



* scos2x + Is 2 cos4x\ e 2l{k - p - 23)x 



exp [ — —s cos 2a; 



p'h- P -2j (~) + lp(Ik-p- 



s \ / s 

2j+l I _ o + h-p-lj-l 



+ ^(k-p- 2j) (lk-p-2j+l (~|) - Ik-p-2j-l (~| 



Hence, for even fc + n, 



Wfc,n(s)-L P ^r 2^ x 



j=max(0,^-p-Lf J) 



ji(^-P-i)! 



(219) 



' P 2 4-p- 2j ' (-|) + fp (4- P -2,+i (-1) + Jfc-p-aj-i (-§)) N 

X 

\ +f(ft-p-2j)(j fc _ p _ % - +1 (-f)-/ fc _ 1> _ 23 -_ 1 (-f)) y 

Again, this is real for real s, so Tik.n — Ti n .k- In particular, for k — n we obtain 



.„ 2 min(LfJ,w-p) 



E 



(n-3-l)\(p + j-l)\ 



sPp\ j\{n-p-j)\ 

j=max[(),n-p- J 



(220) 



/ p 2 /„-p-2,- (-|) + fp (/n-p-2j+l (-|) + 4-p-2i-l (-|)) \ 



(221) 



\ +1 (n-p-2j) (i n _ p _2j+i (-§) - i„- P -2j-i (-§)) y 

We leave odd fc + n as an exercise for the reader. 

Ultimately, Tik.n =/= for k =/= n can be traced back to its origins in the inhomogcncities in the dual 
eigenvalue equations, as was also true for A4,n- A similar calculation for the metric, instead of the dual 
metric, would give 



H k ,n {s) = — T, R(x,p;s)H*f k . n (x,p) = k 2 M k , n (s) = n 2 M n (s) 6 k , n 

27T Jx,p 



(222) 



at least in those cases where is it is well-defined. Hence, for definite energy WFs, we would obtain the 
expected averages: (H) = Ti n (s) /Af n (s) — n 2 . Also, for a linear combination of WFs corresponding to 
the wave function ip = aJ k + (3J n , as in (|207|) above, we would have the more intuitive result: 



w («) = ^ i R fop; *)H*f fop) = \a\ 2 k 2 N k {t 



\ 2 n 2 N n {s) 



(223) 
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We have again assumed that neither k nor n vanishes. 

All this raises the question of how physics can be extracted free from any metric dependence, in general, 
especially in situations where calculations are performed using dual functions that obey inhomogeneous 
equations. Unfortunately, at this time we have no model-independent answer to this question. (However, 
for a brief discussion of possible (semi)classical behavior in the specific exp (2ix) model, please see the last 
subsection in Section IV of [8].) 

12 Conclusions 

There are several types of functions involved in the analysis of imaginary Liouville quantum mechanics: 
Energy cigenfunctions, their duals, and free-particle plane waves related to the eigenfunctions by various 
equivalence maps. By cross-breeding pairs of these functions we have produced a wide variety of hybrid 
Wigner transforms that exhibit a rich diversity of star product relations. Through the use of sums of 
these WFs, and star products of those sums, we have also constructed various phase space equivalence 
transformations and metrics. We have further shown how these phase space kernels can be built from direct, 
elementary solutions of their controlling star product equations. Finally, we have discussed, albeit briefly, the 
import of the metrics for the calculation of physical expectation values. Later, we hope to extend the analysis 
of this paper to include the general class of Hamiltonians of the form H = (p + v) 2 + J2k>o ^ k ex P (ikx). 
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Appendix A. The free particle limit 

This is achieved by taking m — > 0, but first we have to rescale the Bessels and the Neumann polynomials to 
have sensible limits as m vanishes. So we take the biorthonormal set to be {m n A n (me") , m~ n J n (me lx ) } 
which becomes {2 n nle~ mx , 2^\ emx } as m ~ * 0- Indeed, the 2 n n\ factors look weird, but we must attribute 
that convention to Bessel and Neumann. 

Now, what happens as m — > to the metric on the space of dual functions as in [5]? After rescaling the 
Bessels, it is just 

J( X) y) = ]Tm-"J„ (me- fa ) m*J n (me*) ^ £ ^^""W^" = J ° O^"^) (224) 

n=0 n=0 

Another Bessel function. This integral kernel is not unity (i.e. a Dirac delta) but a projector. 

When integrated with (dual) functions that are constants or negative powers of e ly , and their conjugates, 
it gives the expected orthonormality. 

j [ dxdy m n A* n (me vx ) J (x, y) m k A k (me iy ) = 5 n>k (225) 

(2ir) J 

even as m — > 0. But when integrated with positive powers of e ly (not constants) it gives zero. 



Appendix B. Arbitrary functions acting on WFs through star prod- 
ucts 

For any two functions tp and </>, we define the Wigner transform somewhat unconventionally (we do not 
conjugate </>) as 



(x, p) = _ / V (x - y) e 2i «> (x + y) (226) 



where iV will be chosen for convenience, later. Then for any other function G on the phase-space, we have 
the star product 

G (x, p) * (x, P )=g(x,p- ^id^j ^- J e-^ye 2 ^ (x - y) <$> (x + y) dy 

= ^ / e2lVPG ( x -y>P- ^ (x-y)<t>(x + V) dy (227) 

Similarly 

fa (*, p) * G (x, p) = G (x, P + ^)^/ e^e 2 ^ ( x -y)cj>(x + y) e y dy 

f e 2 ^G (x + y,p+Ut\i){x-y)<t>{x + y) dy (228) 



2 N ir I ~ " V~ ' ' 2 



4.3 



For two classes of special cases things simplify considerably here. 
In the first class, if G(x,p) = V (x) then 

V (x) * Ut> ( x >P) = 2^fv(x-y) 4' (x - y) e 2i ™<f> (x + y) dy (229) 
fa (x,p)*V(x) = ^[i>(x-y) e 2i ™V (x + y)<t>(x + y) dy (230) 

Since V (x) here is arbitrary, there is a simple and direct application of (|229p and (|230p to write Wigner 
transforms themselves as star products of wave functions and their Fourier transforms, with an intercalated 
exponential. Namely 

(x,p) =tp(x)* /i0 (x,p) = 4) (x) * -^j- J e 2lyp (j) (x + y) dy 
= ^ I i'(x)*(e- 2ix ^(-2p)) 

= ^- I ^(x)*e- 2i *v*j(-p) (231) 



fa (x, p) = fa (x,p)*<j> (x) = J xjj (x - y) e 2lyp dy * <j> (x) 
-(e 2lxp ${2p))*<j>{x) 



2 N- 
1 



->N- 



T $(p)*e 2ixp *ct>(x) (232) 



where the momentum space wave functions are defined by 

$(2p)=^- J 4> H e~ 2lwp dw , 4> (-2p) = i- j 4>{w) e 2lwp dw (233) 

Reasoning along these same lines also leads to a compact expression for Wigner transforms in terms of star 
products of wave functions with an intercalated momentum-space delta [5], either Dirac or Kronecker. 

fa (x,p) = 2 1 - W V (x) * S 2p * <\> (a:) (234) 

where 5 2p = 6 (2p) = ± e 2lsp ds for x G K and (5 2p = (5 2 p i0 = i e 2isp ds for ifS 1 . 

In the second class, if G(x,p) = K (p) then integration by parts assuming no surface contributions 
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(IPANS) gives 

K (p) * (x, P) = ^j z 2lvP K (p -±id^Tp(x-y)<f>(x + y) dy 



1 



2 N ir 
1 



(x + y) K (~i (K-K)^(x~y) dy (235) 



fa ( x ,p)*K(p) = ^ r -J e^yPK (p +±id^j4,(x-y)0(x + y) dy 

,?ivv K \-~idy + \idA ip(x-y)<f>(x + y) dy 



2 N ir 

f e 2 ^ (x-y)K Q» + flj) J (z + y) (236) 



2^71 

More generally, some care is needed to unravel the remaining derivatives in (|227p and (J225J). For general 
G(x,p) there will appear an operator matrix element involving the operator obtained from G(x,p) by the 
Weyl correspondence. 

As an illustration, consider G (x,p) — xp. Then 

xp * fa (x,p) = -^j- J e 2iyp (x -y)[p- \ id ^J $ ( x -y)<P( x + y) d V 

^ J (x- y) (-\id v (e 2w ) - \e 2 ^id}j iP(x-y)^(x + y) dy 
(~i + (x - y) ii (d y - d x jj tp (x-y)<f>(x + y) dy 



2 JV„ 



upon IPANS. We recognize part of this as the expected operator matrix element [x — y) |i (d y — d x ) ip (x — y) 



2' 

(x — y \ xp \ tp). The other term in the integrand, is the price we must pay because on the LHS we took 
just the simple product, xp, and not the star product, x *p = xp + ^i. Had we taken the later instead, we 
would have 



fa (x,p) = J e 2iyp f(x-y)^i (j^-lt^-ip(x--y)<f)(x + y)dy 

e 2iyp (x-y\xp\tp) <f>{x + y) dy (237) 



X-kp-k 

Z 

Actually, this illustrates the general situation. If the WF is acted on by a "star function" G* (x,p) of x 
and p, by which terminology we mean G* (x,p) is a sum of multinomials of star products of x and p, then 
the result involves an operator matrix clement involving an operator function G (x, p) with the operators 
arranged in the same order as the various x and p factors were arranged in the multinomials of the star 
function. Indeed, that operator is just the Weyl correspondent of the star function. 

g* (x, P ) * jw (*,p) = 2^J e2iyp ( x - y\ G ( x ' p) M ^( x + y) d v ( 238 ) 
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A similar simple operator correspondence applies to acting on the right with G* (x,p), provided we complex 
conjugate the second function. 

(x, p) * G* (x, p) = ^- J e 2 ^ (x - y) (0| G (x, p) \x + y) dy (239) 

Appendix C. Conventions and star product compositions of bilinear 
Wigner transforms 

For functions ip and 4> either on R or on S 1 , we define the Wigner transform as 



( x ,p)= I ^(x~y)^(x + y) e 2 ^dy (240) 



f 

2^7T 

where N is chosen to simplify the normalization of two such functions: TV = for R; N = 1 for S 1 . With 
these choices 

{f +0 ° fibcb (x, p) dp for ifR 
J "°° #l j (241) 

Note that for position coordinates on the circle, x,y € [0,27r], the momentum sum is over all semi-integer 
p so that the complete periodic Dirac delta is produced, X^pez e 2lyp = 2ir5 (y) = 2ir5 (y + 2ir) , such that 
Jo* $ (y) dy = J 2 ^ e S (y) dy = 1 for all e. This is important since in principle the contributing Fourier mode 
numbers of tjj and 4> could differ by cither even or odd integers. 

With these choices for c we also have a uniform appearance for the * composition of two such functions 

U<j> (x, P) * f vx ( x i P) = (0j v) Ux ( x ' P) ( 242 ) 
for both cases. On the RHS we have used the notation 

(<M) = ^/ cj>(x) V (x)dx (243) 

^ JR or §1 

That is to say this metric is local, K (x, y) = ^-5 (x — y) in the language of biorthogonal systems as used in 
the main text, Eqn (j2"0]) . 

For R the result (|242|1 is obtained as in the Overview of [23]. For S 1 the result is obtained with a 
modicum of novelty as follows. Applying the star product directly to each of the integral representations 



4G 



for ftjjrj, and f vx , and carrying out the requisite variable shifts on the integrands, we obtain 

If o.,.„ 1 



U<p {x,P)*fr,x ( X >P) 



1 



2 N ir 



il> (x - yi) <f> (x + j/i) e 2 ^Pe W1 ^ 



2^ 



(2 n tt) 



ij)(x — yi- y 2 ) <P (x + yi- V2) V (x + yi - y 2 ) X ( x + Vi + 2/2) e 2 ^ Vl+V2)p dyidy 2 



(2 N w) 



^ ip (x - yi - y%) x (x + yi + 2/2) e 2l(vi+V2)p d (yi + y 2 ) ^ / <f> {x + y x - y 2 ) rj (x + yi - 2/2) d (j/i - y 2 ) 



(244) 



Now we only need worry about the regions for the two integrations. When the original y± and y 2 are 
integrated over the whole real line, then so are y\ + y 2 and y\ — y 2 , and the result (|242p follows. When 
the original y\ and y 2 are integrated over the circle, it may seem that the y\ — y 2 region of integration 
depends on y\ + y 2 . However, this is actually not the case for periodic integrands. The original coupled 
integrals dy\ dy 2 ■ ■ ■ may be split into uncoupled integrals d (j/i + y 2 ) • • • x ^ ^ — 2/ 2 ) ' ' ' 
upon making use of this periodicity. 
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gions (A, B, and C in the Figure). For 27r-periodic 
integrands the result of the double integration is un- 
changed if the left-most triangle (A) is translated 2ir 
to the right (to become triangle E), and the lower 
triangle (B) is translated up by 2tt (to become tri- 
E \ angle D), to form a 2tt x 47r rectangle (C+D+E) 

on the (yi + y 2 ) x (yi — y 2 ) integration planer 3 
Then periodicity is again exploited to evaluate 
^jll^d(yi~y 2 )(l)(x + yi-y 2 ) r)(x + y 1 -y 2 ) = 
The original 2n x 2tt square representing the region 2tt (0, 77) and hence to obtain the final result (|242j) . 
of integration on the y t x y 2 plane may be partitioned This result is used extensively in the next Appendix, 
by diagonal segments into three right-triangular re- 





C 


A 






B \ 



Appendix D. Non-diagonal WFs, conventional and hybrid 

Consider any biorthogonal system involving countable pairs of functions and their duals, {^mXn}, which 
has an equivalent hermitian system involving corresponding functions, {</>„,(/)„}, orthonormal in the usual 

13 The scale of the final rectangle in the Figure needs to be increased by V2 to correctly display the relative area of the initial and 
final regions on the y\ X j/2 and (j/i + 3/2) x (yi — V2) planes, as is evident from the Jacobian \d (3/1 + 3/2, yi — 2/2) jd (yi , J/2 ) | = 2. 
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sense. That is, in the notation of the previous Appendix, 

(Xk, 4>n) = h,n = (0fe, <t>n) (245) 

In phase space we define the following conventional and hybrid Wigner transforms. 

1 r 271 

e k . n (x,p) = — I (f> k (x- y) 4>n {x + y) e 2iy Uy (246) 



(i 



/*,„ (s,p) = i- / Vfc {x - y) (x + y) e 2lyp dy (247) 



(i 



f k>n (a;,p) = i- / xm (x - y) Xn (x + y) e 2im dy (248) 



5fc ,„ (x,p) = — / ip k (x-y) Xn (x + y) e 2iyp dy (249) 
27r Jo 

h k . n ( x ,p) = JL [ i, k (x - y) K (x + y) e 2lyp dy (250) 

Z7T 



(I 



1 



2jt 



h k , n (x, p) = t- fa (x - y) X n (x + y) e 2lyp dy (251) 

Z7T 







Clearly, in the notation of the main text, e„ (x,p) = e„.„ (x,p), etc. 

From the result in the previous Appendix, (|242| . mixed products are also given immediately by (|242D . 



we immediately obtain the following star products. 


gj.k * hi jTn 


— hj.m $k.l 


(259) 






hj.k * &l,m 


= hj >m 8k,l 


(260) 






Gj.k ^ h>l.m 


— faj.m &k.l 


(261) 


^j.k * &l,m — ^j,7n ^k,l 


(252) 


foj 7 k * h>l.m 


= 9j,m &k.l 


(262) 


fj,k * fl,m = 9j,m $k,l 


(253) 


fj 7 k ^ h>m.l 


foj.m &k.l 


(263) 


fj.k * fl,m = 9m,j $k,l 


(254) 


hk,j * h>i im 


= fj.m 8k, I 


(264) 


fj,k * 9m, I = fj.rn $k,l 


(255) 


hj.k * fl.rn 


= h m ,j 8k,i 


(265) 


9j,k * fl,m = fj.rn Sk,l 


(256) 


hj.k * gi.rn 


= hj >m Sk,l 


(266) 


9j,k * 9l,m = gj.m §k,l 


(257) 


fojik ^ h>l.m 




(267) 


gi.m * fk.j = fm.j Sk.l 


(258) 


There are a few other obvious relations 


that fol- 



low from complex conjugations, using e k . n = £n.k, 
fk.n = f n ,k, and f k .n = fn.k, but in writing these it 
None of these mix the biorthogonal {tpk,Xn} sys- should be kept in mind that the hybrid conjugates, 
tern with the hermitian {</>£,</>„} system. But those g, h, and h, are in general different functions, inde- 



pendent of <?, h, and h. 

Moreover, there are other star products, such as 
3j,/c*/i,m or fj,k*gi,m, which are not proportional to 
5f~ i as simple consequences of the biorthonormality 
of {tpk, Xn} or of {<fik,<j>n}- These involve the scalar 
products {xk,xi), (4>k,4>i), (<f>k,i>n), and (<t>k,Xn), 
which are model dependent and not single Kronecker 
deltas. Thus their evaluation depends on specific 
details of the system. Additional model dependent 
structure exists, such as the form of specific ipa as infi- 
nite series of %s, and vice versa, or the form of specific 
</>s as infinite scries of "0s, and vice versa, which per- 



mit re-expressing some star products as infinite sums 
of others. These re-expressions may amount to rela- 
tions already known to hold in general, as above, but 
often they are new results and highly model depen- 
dent. 

In summary, the full algebra of Wigncr transform 
star products for a biorthogonal system is quite com- 
plicated and varies from one model to another. But 
the subset of products that we have explicitly shown 
are valid for any combined \{4>k, Xn} , {^>fc, <^n} } sys- 
tem. 



Appendix E. A brief overview of biorthogonal systems and density 
operators 

As noted in the introduction, for general biorthogonal systems the Hilbert space structure is at first sight 
very different than that for hermitian Hamiltonian systems inasmuch as the dual wave functions are usually 
not just the complex conjugates of the wave functions. Still, we may keep most of the compact Dirac 
notation for a state, \ip), and its dual, (ip\, provided we just allow (tp\ x) = (x \ip) ^ (x \ip). 

As an alternative notation, we will often write for the dual wave function (ip\ x) = x (x) (no overbar on 
xl) thereby reducing the clutter of symbol decorations. So then, orthonormality for a discretely indexed 
biorthogonal set of states and their duals |(^ r i|| (i- c - a countable basis) reads (ip n \ ipk) = $n,k 

j dx Xn (x) Ipk (x). 

Let us now go through some purely formal manip- wave functions we can always "double-up" the states 
ulations to get a quick overview of the abstract con- and their duals, 
structions that are possible for biorthogonal systems, write the obvious relation \ip n ) = IV'fe) ^fc,« = 
If we are dealing with a discrete biorthonormal basis, 
||"0fc) j ("0n||, then through the usual conjugation of 

14 Kecp in mind that the total number of independent Hilbert space states (L 2 functions) may not have actually changed by 
this doubling-up procedure, e.g. as is certainly the case if the original wave functions and their duals were real functions. But 
this doesn't affect the formal constructions to follow. 



l\^n)A^k\\ 



to include < \ib„), (f/'fcl > and hence 
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S/c I'^fcXV'fcl V'n}- That is to say Now, wc emphasize that it is quite possible, and of- 

ten the case for non-hermitian systems, that R =/= 



\fl> n ) = R\1>n) (268) 

where the operator R is just the formal hermitian 
sum 

R = Yl Vh)^h\ = Rt ( 269 ) 



1. Moreover, it is also often the case that Sk. n ^ 



(ipk\i>n) = Jdx \k (x) Xn (x) and S k>n ^ {ip k \ip n ) 



J dx tp k (x)ip n (x). 

For emphasis and clarity, let us restate the previ- 



This is a positive operator, certainly, but without fur- ous n i ne relations in terms of wave functions ip (x) = 
thcr assumptions about the completeness of the \ip k ), (x their dual functions (ip\ x) = x (x), and the 



R may well annihilate some states, so it is not nec- integral kernel in the position basis 
essarily positive definite. We also have the equally 
obvious relation 

We have 



R(x,y) = (x\R\y) (277) 



(1> n \ = (il> n \R (270) 
Thus the biorthonormahty of the system may also be \n ( x ) ~ J R{x, y) ipn (j/) dy (278) 



written as 

6 n ,k = {i>n\R\1>k) (271) 

The operator R plays the role of a metric in the 
span of {iV'n)}- Given the requisite convergence 
of the sums, for any linear combination of the basis 
states 

k 



R(x,y)=J2xk(x)xk(y) = R(y,x) (279) 



Xn(x) = / i> n (y)R(y,x)dy (280) 



S n ,k = I tpn{x)R{x,y)il>k{y)dxdy (281) 
i>{x) =^c k iP k (x) (282) 



we have a manifestly positive definite norm ^ ^ = J ^ x) ^ {x) ^ (2g3) 

5> fc | 2 = (<A|i?hA) (273) r 

k = / x(xH(x)dx (284) 



Equivalently, we may write 

Icfcl 2 = (7| i>) (274) X {x)= f J{y)R (y, x) dy 

k 



where the dual corresponding to \tp) in the given basis — c„x„ (x) (285) 

n 

is just 



(il>\ = (i ) \R=Y J ^n\ (275) 

n 

as well as the equation dual to this 



and finally 



X( x ) = j R(x,y)ijj(y)dy 
hP)=RW')=J2cM (276) =Y J ^^Jx) (286) 
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H\E)=E\E) , (E\H = E(E\ + (I E \ (293) 



The role of the states and their duals is inter- illustrate more general features of non-hermitian sys- 
changed upon constructing the dual metric terns, so we will focus on them in this discussion. 

~ i / \ / / i ?;+ /r,o^\ Perhaps the most interesting feature of these cases 

R = 2^ lw (w = R ( 287 ) 

k is that the wave functions dual to energy eigenfunc- 

In terms of this operator we have 

tions do not also obey homogeneous equations when 
(VI = (i'\ R = ^"1 ( 288 ) acted upon by the Hamiltonian. Rather, in Dirac 

n 

, , . ~. — r v-^ , notation, the structure is like this: 

|V> = #|V) =]Tc fe |V fe } (289) 

k 

so R is effectively the inverse to R, and vice versa, at 

least on the appropriate subspaces. Correspondingly, The inhomogeneity (I E \ varies with E, usually, but 

we have the formal relations f° r the situation of interest to us here, it is orthogonal 

~ ^-^ , , ; — , , , to all the energy eigenstates, (Ie \E') = 0. For the 

RR = J2\<Pk)(i>k\ (290) 6 ' Ul ' 

ft PT symmetric theories of interest, the energy eigen- 

RR = ^\i>k){i>k\ (291) values are real. 

It is evident at this point that we may have some 

freedom in our construction of density operators, 

The first of these acts as the unit operator on the span , . , r~r\ , , , , mi • • 

given that we may nave \%pk) \ipk)> etc. Inis is 

of II il'k)} while the second acts as the unit on the span „. „ , . . , „ 

u so. I here are lour choices m general. _bor pure 

of I \il>k) <■■ So RR and RR are at least projection . ,. . rj.ii.- 

L 1 'J J states m the various spans ol the basis vectors, as 

operators onto those respective spaces, but they are .„ , , , . , .. , 

specified above, but otherwise arbitrary, we have 

not necessarily 1 acting on all functions. The effect 

f 1 P=IV)(VI=P t (294) 

of RR acting on the span of s IV'fe) f i s n °t clearly 

discernible, nor is RR acting on the span of {iV'fc)}- ^ _ I^XV 1 ! — P (295) 

For many PT symmetric systems, the duals are as we u as the less symmetrical 

proportional to the wave functions. When these arc — . 

(~p) = |V) (VI = (Pl f (296) 

chosen to be PT eigenstates, as is usually the case, 

(p~) = |V) (V| = (» f (297) 

then in a broad class of situations Xn (%) = iVn (^)- 

Again, we stress the absence of complex conjugation. These are all interrelated by the R and R operators. 
However, there are some interesting situations where / ~\ _ ^ _ (298) 

the dual functions are not so simply related to the 



fe 

RR = RRRR , RR = RRRR (292) 



wave functions [51 [5] . This is particularly the case 
when dealing with eigenfunctions at so-called "spec- 
tral singularities." These exceptional cases seem to 



(» = Rp = pR (299) 
p = RpR={p~)R = R(~p) (300) 
p= RpR= (~ p)R = R{p~) (301) 
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When (-01 0) = 1, the various density operators also 
obey several variants of the standard pure state con- 
dition. 



PP = 


oo 


(302) 


pp = 


(» 


(303) 


0000 = 


oo 


(304) 


(»(» = 


(» 


(305) 


Pip') = 


p 


(306) 


(~p)p= 


p 


(307) 


(p~)p = 


p 


(308) 


p(~p) = 


/> 


(309) 



However, we note that p 2 ^ p and p 2 ^ p, in general, 
but rather 

pRp = p , pRp = p (310) 

as implicitly stated in the previous set of equations. 

The results of the previous Appendix can be re- 
cast into the density operator language. A full 
equivalence with the results there can be achieved 
provided there is another set of states and their du- 
als, {\4>k) ■ (071 1 }, in one-to-one correspondence with 
jj0/c) , (0 n | |, where the other set is equipped with 
a trivial metric (i.e. R = 1), as would be the case 
for a biorthogonal system which admits a similarity 
transformation to an hcrmitian Hamiltonian. The 



full equivalence is then given by 

efe,„ ~ |0fe) (0„| (311) 
/fc,n~hk)(lM (312) 

JZn~Vh)hK\ (313) 

g k ,n ~ l^fe) Wn\ (314) 

hk.ri ~ U'k) (0n| (315) 

Kn - \<t>k) (0^1 (316) 

The star products of the previous Appendix corre- 
spond to the operator products of these dyadics, in 
an obvious way. 

When the metric R is positive definite and hence 
invertible, the other set of states could be realized 
through its hcrmitian square root S, with 

R = S 2 (317) 

In this case we may take 

|0 fc ) = S U> k ) (318) 

(0^| = (0„|5 (319) 

100 = S- 1 |0 fe ) (320) 

(0„| = (SKIS- 1 (321) 

If R is positive but not definite, only the first two of 
these are guaranteed valid, but then the role of S" -1 to 
give the latter two relations may be effectively played 
on the appropriate subspace by the "other" square 
root, S, where 

R = S 2 (322) 
\il>k) = S\4> k ) (323) 
(0„| = (0J5 (324) 
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Appendix F. Operator expressions from Weyl transforms 

Real Liouville theory (adapted from [Til]) Given the factorized phase-space generating function 

g(z;x,p) = K lp (e z ) exp (-1 e 2 ^ (325) 
what is the operator corresponding to it? According to Weyl's prescription the associated operatoil 15 ! is 

G(z) — [ drdadxdp Q(z;x,p) cxp(ir(p — p) + ia(x — xj) (326) 

(27T) 2 J 

= - — / drdadxdp exp(irp + iax) exp ( e 2x ~ z — iax) Ki p (e z ) exp(— irp). 

(27r) J \ 2 J 

The integrals over x and p may be evaluated separately, if the a contour is first shifted slightly above the 
real axis, a — ► a + ie, thereby suppressing contributions to the cc-integral as x — ► — oo. Now s = \e 2x ~ 
gives 

J + °° dx exp (~ e 2x ~ z -i(a + ie) x\ = ^ {2se z )- l{,J+le)/2 exp (s) 

= - e - liz+ln 2)ff / 2 r (-i (a + ie) /2) (327) 

But then again 

dp K lp (e z ) cxp(-irp) = - / rfX e' e ' coshx 2ttJ (X - t) = tt e" e * coshr (328) 

-OO ^ J — OQ 



So 



G(z) = / drda e -*(*+ln2) CT /2 r (-i (<r + ie) /2) e" e " cosllT cxp(irp + iax) (329) 

87T J 



The shifted er contour avoids the pole in T at the origin. 

Re-ordering with all ps to the right (thereby departing from Weyl ordering but without actually changing 
the operator) yields exp(irp + iax) = cxp(icrx) exp (iar/2) exp(irp). Performing the a integration before 
the t integration permits taking the limit e — > to obtain 

G(z) = i J dT (^J daT(-i(a + ie) /2) cxp(iax + iar/2 - ia (z + In 2) /2)J e" e * cosllT cxp(irp) 

= -L /" dT (47TCXP ^_ e 2x+r-( Z +ln2)^ g-e* coshr exp ( irp ) 

i / dre- eZcoshT exp f -ie 2x+r - z ) exp(irp) (330) 



2 J V 2 

This is the operator correspondent of Q(z; x,p). In this form it is straightforward to show that 

[p 2 + e 2x ,G(z)] =0 (331) 



3 The operator corresponding to O will be distinguished from it by using the same letter but sans serif: O. 
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for all values of z, which suggests an interpretation of G(z) as a propagator and a concomitant interpretation 
of z as a generalized "time." 

The form (|330p also leads to a more intuitive Hilbert space representation. Acting to the right of a 
position eigen-bra, we have (x\ x = (a;| x, while the subsequent exponential of the momentum operator just 
translates, (x\ exp(irp) = (x + r|. So the full right-operation of G is 

(x\ G(z) = l -jdr{x + t\ cxp (-\e 2 ^-* _ ^ _ ^-r^ 

= \Jdy (y\ exp (-\e x+y - z - \<?^*°) (332) 

Inserting 1 = J dx \x) (x\ gives G(z) = J dx \x) (x\ G(z), and leads to a coordinate space realization of the 
operator involving an x, y-symmetric kernel, in which form it is clear that G(z) = G(zY for real z. 



G(z) = \J dxdy \x) exp {—f^" ~ \e x ~ 



(y\ (333) 



All the operator characteristics are now carried by the dyadics \x) ■ ■ ■ (y\. The composition law of this 
operator parallels its phase-space isomorph. 

G(u)G(v) = ^Jdw exp ^-i (e u+v - w + e u - v+w + e - u+v+w )^ G{w) = G(v)G(u) (334) 

Imaginary Liouville theory Let us apply these same methods to the imaginary Liouvillc case. Here, 
we can just write down the final answer from what we know about real Liouville theory, as given above, 
and then work backwards, a procedure often followed in [S]. It is easily seen that the following operator 
similarity transformation is hermitian, for real s, and converts operator H into operator H^. 



F(s) = - ( dxdy \x) cxp ( — e ^- x+vS > - -e 1 ^^ - '- e ^+v) ) ( y | (335a) 
— I V _^ -S' 2 2 



F(s)H = H t F(s) (335b) 



H = p 2 + e 2ix , H t = p 2 + e 2lx where x, p are operators (335c) 

However, it is not obvious that F (s) is positive definite, even on a restricted range of p eigenstates. 

Acting on a position eigenbra, (x\, it can be seen that F is given by the parity operator multiplying 
a single parametric integral, whose integrand can be written in factorized form as a product of a spatial 
translation and an exponential of a momentum translation, just like the real Liouville case. That is to say 

F( s ) = 1 P J da e - scos(T cxp ^i-e^+^J exp (tap) (336) 
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where again x and p are position and momentum operators, and where P is the usual parity operator: 
(x\ P = (— x\, PxP = —x, PpP = —p. The range of the a integration is whatever is needed for completeness of 
position eigenstatcs: 1 = J da \a) (<r|. The form (|336[) should be compared to (|330p . with the identification 

,s = — e z . 

The presence of P in F(s) has interesting consequences. From (|336p it follows for real s that we can write 

F (s) = J da e - scosa D (s) P D f (s) (337) 

where 

D( S ) = cxp(-^ap)cxp(^e- 2ix ) , D t (s) = exp Q-/^ exp (\iapj (338) 

From the form (|337p it would seem that F is not positive, although it is manifestly hermitian. From the 
form p36|) we can also check that F = for real s, but not so easily as from <|335a[) or from ()337|) . 

In the course of obtaining (|336[) from (|335ap it is remarkably easy to see how the parity operator appears. 
No prior knowledge of its presence is needed. However, if we remove P from the RHS of (|336p . we see that 

[H,PF(.s)] = = [H t ,F(.s)P] (339) 

which is similar to (|331|) . All this again suggests an interpretation of PF(s), or F (s)P, as a propagator for 
the system governed by H, or H^, as well as an interpretation of s as a generalized time. Moreover, from 
(|336p the phase-space correspondent of PF (s) is obviously given by 

T{s-x,p) = \J d ° e~ scos,T exp ^l-e l( - 2x+a ' ) ^ * exp {tap) (340) 
= iy dae- SCOBa cxp(iap) exp (^-e 2lx ^j (341) 

This is nicely factorized into a function of p times a function of x. For integer p, with a integrated over 
[0, 2tt], it is simply 

F{s; i,peN)= nip (-s) exp (j-/ 1 *) ( 342 ) 

This should be compared to (|325|) . again with the identification ,s = — e z . 

We leave as exercises for the interested reader to determine the operators that result from taking the 
Weyl transforms of the phase-space metrics and/or the dual metrics. 

FTeX 
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